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PART  ONE 
GENERAL  THEORY 


I.  INTRODUCTION 

The  problem  of  electromagnetic  excitation  of  a  wire  through  an 
aperture-perforated  conducting  plane  has  been  studied  in  several  papers. 
Kajfez  [1]  derived  the  equivalent  sources  and  traveling  current  on  the 
wire  excited  by  a  small  circular  aperture,  but  he  neglected  the  interaction 
of  the  wire  on  the  electric  field  in  the  aperture.  The  axial  distribution 
of  magnetic  current  on  the  slot  is  sometimes  markedly  affected  by  the 
presence  of  the  wire.  Butler  and  Umashankar  [2]  considered  the  boundary- 
value  problem  of  thin  wire  of  finite  length  behind  a  slotted  conducting 
palne,  obtaining  a  coupled  set  of  integro-dif ferential  equations,  which 
they  solved  numerically.  A  later  paper  [3]  by  Umashankar  and  Wait  treated 
an  infinite  cable  placed  behind  a  slot-perforated  screen  problem  by  a 
Fourier  transform  method.  However,  there  are  errors  in  their  equations 
and  in  their  computed  results.  For  example,  the  term  cos  Az  is  missing 
from  their  Eq.  (19d),  and  no  transmission  line  component  of  current  is 
apparent  in  their  Fig.  5. 

In  this  report  we  obtain  the  functional  equations  for  the  problem 
using  the  equivalence  principle  [4],  and  reduce  these  equations  to  matrix 
form  by  means  of  the  method  of  moments  [5].  The  matrices  are  recognized 
as  generalized  network  parameters,  such  as  voltages,  currents,  admittances, 
and  impedances,  as  described  in  [5]  and  [6]. 

We  next  specialize  the  equations  to  the  narrow  slot  and  infinite 
thin  wire  problem.  As  electric  current  expansion  functions,  we  take 


2 


outgoing  transverse  electromagnetic  (TEM)  currents  traveling  on  both 
semi-infinite  halves  of  the  wire,  and  triangle  functions  to  represent 
the  higher  order  modes  in  a  finite  region  near  the  aperture.  This 
results  in  a  finite  sized  matrix  for  the  infinite  wire,  so  that  the 
solution  is  similar  to  that  for  a  finite  wire. 

Finally,  in  addition  to  excitation  by  a  plane  wave,  we  con¬ 
sider  excitation  by  a  TEM  traveling  wave  on  the  wire.  The  perturbation 
due  to  the  effect  of  the  slot  can  be  obtained  from  the  outgoing  TEM  cur¬ 
rents  in  both  directions.  The  equivalent  circuit  parameters  for  the  TEM 
wave  are  then  calculated  from  these  outgoing  traveling  waves.  Once  the 
equivalent  circuit  parameters  are  obtained,  we  can  calculate  the  currents 
traveling  on  the  wire  for  the  case  of  a  wire  terminated  by  arbitrary 
loads. 

II.  FORMULATION  OF  THE  PROBLEM 

Figure  1  shows  the  problem  to  be  considered  and  the  coordinates  and 
notation  to  be  used.  The  infinite  conducting  plane,  which  has  zero  thick¬ 
ness,  covers  the  entire,  z=0  plane  except  for  the  aperture.  The  aperture 
is  rectangular  in  shape  with  length  b  and  narrow  width  w  in  the  y  and  x 
directions,  respectively.  The  plane  divides  the  space  into  two  regions, 
called  region  a  and  region  b.  There  is  an  infinitely  long  wire  whose 

axis  is  along  the  x-direction  and  at  a  distance  z=d  and  y=y  in  region  b. 

c 

The  excitation  of  the  slot  is  a  uniform  plane  wave  incident  from  the 
region  a.  The  problem  is  primarily  that  of  finding  the  tangential  electric 
field  on  the  slot  and  the  current  on  t he  wire.  Secondarily,  it  is  that 
of  finding  the  parameters  of  the  equivalent  network  seen  by  the  TEM 
traveling  wave  on  the  wire. 


Fig.  1.  Thin  infinitely  long  wire  behind  a  narrow  slot  perforated  conducting 
plane. 


We  use  the  equivalence  principle  [A]  to  divide  the  original 
problem  into  two  parts,  as  shown  in  Fig.  2.  The  field  in  region  a 
remains  unchanged  if  the  slot  is  closed  by  a  conductor  and  the 
equivalent  surface  magnetic  current 

M  =  u  x  E 
—  — z  — 

is  placed  over  the  slot,  where  E  is  the  electric  field  in  the  slot 
of  the  original  problem  and  u  is  the  unit  vector  in  z  direction.  The 

“Z 

incident  wave  _E1,  JJ1  must  be  kept  in  region  a.  This  equivalence  is 
shown  in  Figure  2a.  The  field  in  region  b  remains  unchanged  if  the 
slot  is  closed  by  a  conductor  and  the  equivalent  magnetic  current  -M 
is  placed  over  the  slot.  The  electric  current  I  on  the  wire  must  be 
kept  the  same  as  the  original  problem.  This  equivalence  is  shown  in 
Fig.  2b. 


X 


( a )  Region  a 


X 


Fig.  2.  The  equivalent  problem 


The  use  of  M  in  region  a  and  -M  in  region  b  ensures  the  continuity 
of  the  tangential  electric  field  across  the  slot.  In  addition,  we  have 
two  more  boundary  conditions  to  enforce:  (a)  the  tangential  magnetic  field 
must  be  equal  on  each  side  of  the  slot,  and  (b)  the  tangential  electric 
field  must  vanish  on  the  surface  of  the  wire.  These  two  conditions  give  us 
two  equations  from  which  we  can  calculate  the  unknowns  M  and  I. 

The  electric  and  magnetic  fields  in  region  a  are  given  by 

E3  =  Ea(M)  +  E1  (2) 

H3  =  Ha(M)  +  H1  (3) 

where  Ea (M) ,  Ha (M)  are  the  fields  from  M  in  Fig.  2a,  and  E1,  H1  are  inci¬ 
dent  fields  in  region  a.  Note  that  all  the  fields  are  computed  with  the 
slot  shorted.  Similarly,  we  denote  the  electric  and  magnetic  fields  in 
region  b  to  be 


=  Eb(-M)  +  Eb(I) 

(4) 

=  Hb(-M)  +  Hb(I) 

(5) 

where  12  (-M) ,  11b  (-M)  are  the  fields  from  -M  in  Fig.  2b,  and  Eb(I),  Hb(I) 
are  the  fields  from  _I  on  the  wire  in  Fig.  2b.  Again,  all  the  fields  are 
computed  with  the  slot  shorted. 

To  enforce  the  two  boundary  conditions  (a)  and  (b),  we  set 


Ha(M)  +  Hb(M)  -  Hb(I)  -  -  H J 

over  A 

(6) 

e£(m)  -  e£(i)  -  0 

over  B 

(7) 

wher*.  A  and  B  .tote  the  slot  and  the  wire  surfaces,  respectively.  The 
subscr^t-e  t  denotec  the  tangential  component  over  A  and  B.  We  have  used 
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the  linearity  of  the  operators  in  deriving  (6)  and  (7),  by  setting  = 

an<3  E^(-M)  =  -E^(M).  Equations  (6)  and  (7)  are  .rector  equations  for 

determining  the  unknowns  M  and  _I. 

We  now  reduce  (6)  and  (7)  to  matrix  equations  using  the  method  of 

moments  [5].  We  define  a  set  of  expansion  functions  {M  ,  n=l , 2, . . . ,NA} 

— ■ n 

and  express  the  magnetic  current  as 


NA 

m  =  y  V  M 

—  ,  n  — n 

n=l 


(3) 


where  are  coefficients  to  be  determined.  We  define  a  set  of  expansion 

functions  {I  ,  n=l,2, . . . ,NB},  and  express  the  electric  current  as 
— n 


NB 

I  =  y  a  I 

—  ,  n  — r 

n=l 


(9) 


where  are  coefficients  to  be  determined.  We  substitute  (8)  and  (9)  into 


(6),  and  (7),  using  the  linearity  of  the  operator  and  obtain 

NA  NB  .  NB  . 

-  I  V  H*(M  )  -  I  V  1T(M  )  +  l  a  H"(I  )=H* 

S  n— t  — n  L.  n-t  — n  L .  n  -t  ~n  — t 

n=i  n=l  n=i 

over  A,  and 


NA  NB  , 

I  V  E°(M  )  -  T  a  E, (I  )  =  0 

L .  n-t  — n  L.  n-t  — n 

n=l  n=l 


(10) 


(ID 


over 


B.  We  define  testing  functions  {M  ,  m=l, 2, . . . ,NA)  on  A  and 


(i  ,  m»l, 2, . . . ,NB}  on  B,  and  the  symmetric  product 
— m 


<F,  G>  = 
—  —  s 


G  ds 


(12) 


where  S  represents  the  surface  of  A  or  B.  We  take  the  symmetric  product 


of  (10) 

with  each  M  on  A  to  obtain 
-m 

[Ya  4-  Yb]V  +  [TJI  -  I1 

(13) 

where 

[Va]  -  [<  ^(M,»  A1 

(14) 

(Yb]  -  [<  A) 

(15) 

C'i. 
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The  [Y  ]  and  [Y^]  are  called  generalized  admittance  matrices  for 
region  a  and  b,  respectively.  The  matrix 


[T]  =  [<  M  ,  H  (I  )>  ] 
— m  t  — n  A 


(16) 


is  called  the  coupling  matrix  from  B  to  A.  The  column  matrix 


->i  a  i  _ 

i  =  [<m  .  «i>A] 

—m  — t  A 


(17) 


is  called  excitation  vector.  Finally,  the  column  matrices 


V  =  [vn] 


(18) 


1  =  [an] 


(19) 


are  called  generalized  voltage  and  current  vectors,  respectively. 
Next,  we  take  the  symmetric  product  of  (11)  with  each 


I  to  obtain 
-m 


[T] V  +  [Z]I  =  0 


(20) 


where 


[T]  =  [<I  ,  E  (M  )>  ] 
— m  —t  n  B 


(21) 


is  called  coupling  matrix  from  A  to  B,  and 


[Z]  =  t<-I  ,  E“(I  )>  ] 
-^n  t  n  B 


(22) 


is  called  generalized  impedance  matrix  for  the  wire. 


We  now  have  the  matrix  equations  (13)  and  (20),  from  which 
we  can  solve  for  V  and  I.  Note  that  [Ya]  and  [Y*5]  are  evaluated 
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for  the  slot  as  if  the  wire  were  not  present,  and  are  hence  the  same 
matrices  as  described  in  [6].  [Z]  is  evaluated  for  wire  as  if  the  slot 

were  not  present,  that  is,  with  the  slot  closed  by  a  conductor.  It  can 
be  evaluated  in  terms  of  currents  plus  their  images  radiating  into  free 
space. 


III.  EVALUATION  OF  THE  MATRICES 

Now  we  specialize  the  problem  to  a  narrow  rectangular  slot  and  a 
thin  wire,  and  evaluate  all  the  matrices  in  equations  (13)  and  (20).  We 
use  Galerkin's  method,  which  means  that 


M  =  M 


I  =1 
— m  m 


(m  =  1,2, .. . ,NA) 
(m  =  1,2,. ..,NB) 


(23) 

(24) 


This  will  save  us  some  computations.  In  particular,  from  reciprocity 
[4] ,  we  have 


<M  ,  H. (M  )>  =  <M  ,  H.(M  )>. 

~tn  t  — Ti  A  — n  t  -m  A 


<1  ,  E  (I  )>  =  <1  ,  E  (I  )> 

t  n  B  n  — t  “in  B 


Substituting  (23)  and  (24)  into  (14),  (15)  and  (22),  we  find  that 


Ya 

mn 

«  Ya 
nm 

or 

nra]  -  tfa] 

(25) 

Yb 

mn 

=  Yb 
nm 

or 

[Yb]  -  [Yb] 

(26) 

Z 

mn 

*  Z 

nm 

or 

[Z]  -  [Z] 

(27) 

where  denotes  the  transpose  of  a  matrix.  Equations  (25),  (26)  and  (27) 


indicate  that  the  generalized  admittance  and  impedance  matrices  are 
symmetrical  matrices.  Again,  from  reciprocity,  we  have 
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<1  .  E"(M  )>  = 

t  n  B 


"  ,  H“(I  )> 

t  - m  A 


Substituting  (23)  and  (24)  into  (16)  and  (21),  we  find  that 


T  =  -  T  or  [f ]  =  -  [T]  (28) 

mn  nm 

Hence,  we  only  need  to  evaluate  one  of  the  matrices  [T]  and  [T ] . 

(A)  Evaluation  of  [Z] 

First,  we  have  to  choose  appropriate  expansion  functions  for  the 
current  so  that  our  results  converge  rapidly  to  save  computation.  We 
note  that  the  wire  and  the  complete  conducting  plane  form  a  transmission 
line.  We  therefore  expect  that  the  presence  of  the  slot  on  the  plane  will 
excite  outgoing  traveling  TEM  waves  along  the  wire  plus  evanescent  waves  near 
the  slot.  We  can  expect  that  the  TEM  traveling  waves  will  propagate  to 
infinity  without  attenuation  if  the  materials  are  loss-free.  The 
evanescent  waves,  which  are  higher  order  modes  on  the  transmission  line, 
attenuate  rapidly.  Therefore,  we  assume  that  the  evanescent  waves  will 
exist  only  over  a  short  region  near  the  slot.  Beyond  this  region,  the 
evanescent  waves  are  so  small  that  they  do  not  contribute  significantly  to 
the  field,  and  hence  may  be  taken  to  be  zero.  It  is  therefore  convenient 
to  choose  the  following  expansion  functions  for  electric  current  I  on  the 
wire: 

-oo  <  x  <  oo  (29) 

n  -  2,  3,...  NB  (30) 


Ix  -  e"jk|x|. 


I  -  T  (x) 

*“T\  Tl 
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where  k  is  the  wavenumber  of  incident  wave,  u  is  the  unit  vector  of  x 

—x 

direction,  and  T^x)  are  triangle  functions  as  defined  in  Appendix  B. 

We  can  see  that  L.  represents  the  outgoing  TEM  mode,  and  I  (n  >  1)  repre- 
sent  the  evanescent  currents  near  the  slot.  We  here  assume  that  the  higher 
order  modes  exist  only  in  the  region  x  <  |  L./2 1 . 

Now  we  can  evaluate  the  elements  of  the  generalized  Z  matrix.  As  a 
matter  of  fact,  for  m  4  1,  n  ^  1,  it  is  the  Z  matrix  of  a  wire  whose  length 
is  L  above  a  conducting  plane,  as  is  described  in  Appendix  B.  The  results 
are 


=  j<4iZ 
mn  1 


(31) 


where  Z^  and  Z^  are  expressed  in  (B-18)  and  (B-19) . 
For  m  =  n  =  1,  from  (22)  we  have 


zu  ■  -  'ir  W'. 


(32) 


where  E^(I^)  is  the  tangential  electric  field  on  the  wire  surface  due 


to  outgoing  TEM  currents.  It  is  derived  in  Appendix  A  for  the  field  due 
to  outgoing  TEM  traveling  currents  in  free  space.  In  our  case  there  is  an 
infinite  conducting  plane  near  the  wire.  By  the  image  theory,  it  is  equi¬ 
valent  to  the  total  field  due  to  current  on  the  wire  and  the  current  on 
its  image,  both  in  free  space.  Therefore,  using  (A-14),  we  obtain 


-jkRo  -jkR^ 

Ex*-1*  =  27To>e  (  R  R*  ) 

o  o 

where 

rq  =  V®2  +  x2 

R'  -  fid2  +  x2 

o  f 


(33) 


(34) 


(35) 
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Substituting  (29)  and  (33)  into  (32),  we  find 


-11  27ru)e 


-jkR  -jkR' 


The  integral  can  be  evaluated  as 


-jkR 


00  -jk(R  +x) 
'  e  0 


dx  + 


0  -jk(R  -x) 

*  e  0 

R' 
o 


dx 


°°  -jk(R  +x) 
e  ' 


dx  =  -  2  [Ci(ka)  -  jSi(ka)] 


where 

Ci(x)  = 

Si(x)  = 

Using  (37),  we  can  reduce  (36)  to 


CO 


X 

► 

0 


cos  u 
u 

sin  u 
u 


du 


du 


(36) 


(37) 


(38) 

(39) 


Z11  =  iSc  [Ci(2kd)  -  Ci(ka)  +  jSi(ka)  -  jSi(2kd)]  (40) 

For  m  Mi  n  =  1,  substituting  (30)  and  (33)  into  (22),  we  obtain 


ml  2iru)e 


2m+3 


-jkR 


-jkR’ 


T  00  (S— 

m  K 


R’ 

o 


-)dx 


(41) 


2m-l 


where  the  notation  is  the  same  as  that  used  in  Appendix  B.  We  use  the 
same  approximation  as  that  in  Appendix  B,  that  is,  the  testing  function 
T  (x)  is  approximated  by  four  impulses.  Hence,  (41)  becomes 


ml  2ttwe 


4 

l 

A-l 


C*(A) 

m 


-jkR 


-jkR’ 
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AX 


R’ 

o 


-)dx 


(42) 


2m-2+A 
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Using  Green's  function  ip  as  defined  in  (B-21),  we  can  rewrite  (42)  as 


ml 


=  —  l 

“e  A-l 


C*(A) 

m 


[<K0,  Q 


2m-2+A) 


-  <K0, 


Q ' 2m- 2  +  A* 1 


AX 


2m-2+A 


(43) 


All  the  elements  of  [Z]  matrix  are  now  evaluated  explicitly  in  (31),  (40), 
(43),  and  (27). 

(B)  Evaluation  of  [Ya]  and  [Yb] 

a  K 

According  to  the  definition  of  [Y  ]  and  [Y  ],  Eqs.  (14)  and  (15), 

SL  b 

[Y  ]  and  [Y  ]  are  the  matrices  obtained  by  using  the  expansion  function 
as  a  magnetic  current  radiating  into  region  a  and  b,  respectively,  with 
the  infinite  conducting  plane  present  and  the  slot  closed  and  the  wire  not 
present.  Since  the  two  regions  a  and  b  are  identical  half  spaces,  their 
generalized  admittance  matrices  are  same.  Therefore  we  have 

[Ya]  =  [Yb]  =  [Yhs]  (44) 

lis 

where  [Y  ]  is  the  generalized  admittance  matrix  for  the  slot  opening  into 

tlS 

half  free  space.  Applying  the  image  theorem,  we  find  that  [Y  ]  is  the 
admittance  matrix  obtained  by  using  21^  as  magnetic  current  radiating  into 
free  space  everywhere. 

lis 

Now  we  specify  [Y  ]  for  our  narrow  slot  case.  In  this  problem 
the  width  W  is  small  relative  to  the  wavelength  and  the  length  is  large 
compared  to  its  width.  The  electric  field  in  the  slot  is  principally 
transverse  to  the  longer  slot  dimension  (transverse  to  y),  and  has  a  known 
transverse  variation  given  in  the  literature  [7].  This  transverse  variation 
of  electric  field,  or  its  equivalent  magnetic  current,  is 


J 


13 


f(x) 


=  —  [(-)2  - 

TT  1 


2.-1/2 

x  ] 


(44) 


where  w  Is  the  width  of  the  slot.  Also,  the  magnetic  current  has  only 
a  y  component.  Hence,  the  magnetic  current  expansion  functions  may  be 
expressed  as 

M  (x,y)  ■  f(x)  T  (y)  n  ■  1,2, ...,NA  (45) 


where  vector  functions  T  (y)  are  triangle  functions  defined  in  Appendix  B. 

— n 

Now  we  evaluate  the  magnetic  field  H  (M  )  on  the  slot.  For  H  , 

— t  - n  — t 

we  have  the  following  formulation: 


94> 

H  (M  )  =  -  jcoF  -  ~ 
y  -n  y  9y 


y  4tt 


;JK 


-m 

—  ds' 


4>, 


m  4tt(j  J 


-jkR 

°m  V-  d3’ 


O  =  -  V  •  M 
m  jo)  — n 


(46) 

(47) 

(48) 

(49) 


where  F^  is  the  only  component  of  electric  vector  potential,  and  4>m  is 

the  magnetic  scalar  potential.  The  integral  is  carried  over  the  region 

where  M  exists.  R  is  the  distance  from  the  field  point,  which  is  on 
“II 

the  slot,  to  ds', 

r  =  -\/x'^  +  (y  -  y')^ 


Substituting  this  into  (44)  we  calculate  the  integral  S  as 


T  -Ky-y’) 


f  -  jvfy  -x  v  x  +  (y-y') 


2  2  ~  e-jk  Vx,2+(y-y')2 

T  (y')  -  - —  -  ■  ■  dx'dy* 

I  -  f  V (f)2  -  x'2  Vx’2  +  (y-y’)2 


The  change  of  integration  variable 


i  w  .  a 
:  '  2  sl”  2 


reduces  (51)  to 


-  2  0  a 

,,2  ,  w  .  2  — 
-y  )  +  / —  sin  2 


Tm(y,)dy' 


7r  -jk"Y  (y-y1)  +  4 

e 


J  in2  .2a 

-TT  271'  y  (y~y')  +  X  sin  -r 


,.2  . 

-y’)  +  t- 


-jkV(y-y')  +4 


T  (y’)  ^ -  dy’  *  | 

I  ■  y^TTf-  -  | 


t Jy')  —  dy’  (52) 

b  m  Re 


Here  we  have  used  the  relation 


77  -jkV(y-y’)2  +  4"  sin2  f  -JkV^y-y')2  +  (4? 


1  1/ /  2  w^  .2a 

-tt  V(y-y  )  +  t  sin  T 


r  , .  2  ,w.  2 

(y-y’)  +  (r) 


Re  -  V(y  -  y')2  +  (f)2 


With  this  result,  (47)  through  (49)  can  be  reduced  to 
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F 
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e 

4tt 


b 

2 

r 


t  (y’) 

m 


1 

4  Try 


(55) 


(56) 


a 

m 


p  V*’> 


(57) 


of  (46),  together  with  (55)  through  (57),  is  recognized  to  be  of 

the  form  as  the  magnetic  vector  potential  associated  with  a  thin  wire  of 

radius  w/4.  It  is  the  dual  problem  to  the  formulation  of  [Z]  matrix  for 

a  thin  wire  of  radius  w/4.  Therefore,  using  this  duality,  we  obtain  the 

formula  for  the  [Y]  matrix  of  a  slot  with  triangular  expansion  functions 

T  (y)  to  be 
m 


where 


Yhs  =  2  (jooe  Y,  +  ~  Y_) 
mn  VJ  1  juy  2' 


4  4 

Y1  =  l  l  C^(A)  c£(B)  Ay  Ay  <KQ  ,  Q  ) 
A“1  B-l  1  3  13 

4  4 

Y2  -  l  l  d£(B)  Ay  Ay  <HQ  ,  Q  ) 

A-l  B-l  3  3 


i  -  2m  -  2  +  A 


j  -  2n  -  2  +  B 


The  Green's  function  \p  is  defined  similar  to  (B-21) , 

1 


>KQi»  Qj)  *  4-rrAy 


i 


Pj+1  e"jkRe 


dy' 


rj 


(58) 

(59) 

(60) 

(61) 

(62) 

(63) 


1 


! 


I 

} 

i 

i 

4 


where  is  at  the  center  of  the  ith  segment  Ay^  of  the  slot,  P  is  the 

jth  dividing  point,  and  Rg  is  defined  in  (54).  Note  that  the  coefficient 

2  in  (58)  comes  from  the  fact  that  the  magnetic  field  produced  by  2M  in 

— n 

free  space. 


(C)  Evaluation  of  [T]  and  [T] 

In  order  to  evaluate  [T]  from  (16),  we  have  to  calculate  H^(I  ) 

— t  — n 

on  the  slot.  In  our  coordinate  system,  we  have 


H7 (I  )  =  H  (I  )v 
— t  — n  y  — n  — y 

From  the  geometry  shown  in  Fig.  3,  we  obtain  the  following  relations: 


H  *  u  •  H ,  u 
y  -y  $  -<(> 


u.  *  cos  a  u  +  sin  a  u 
-<p  -y  -z 

where  u  and  u  are  unit  vectors  in  the  y  and  z  directions,  respectively, 
— y  —  z 

u.  is  the  unit  azimuthal  vector  of  the  circle  centered  on  the  wire,  and  a 


Fig.  3.  Wire  and  slot 
in  y-z  plane. 
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Is  the  angle  between  _u^  and  u^.  It  is  given  by 


cos  a  =  d/p 


(64) 


where 


=  t/( y  -  y  )2  +  d2 


(65) 


Therefore 


H  = 

y 


V(y-yc)2  +  d2 


(66) 


For  n*l,  ( 1.^ )  Is  expressed  in  Appendix  A,  (A-13) .  Substitution 

into  (66)  gives 


W  ■ 


-jra. 


^[(y-yc)2  +  d2] 


(67) 


Note  that  (67)  already  includes  the  field  due  to  the  image  of  L^,  and 

the  wire  is  now  taken  in  x  direction  instead  of  in  the  z  direction  as  in 

Appendix  A.  Since  the  axis  of  the  slot  is  in  the  y-z  plane,  we  have 

R  ■  p.  Hence,  (67)  can  be  rewritten  as 
o 


W  = 


-jk  7(y-yc)2-W2 


"[(y-yc)2  +  d2] 


(68) 


Substitution  of  (23),  (68)  into  (16)  gives 


T  ,  =  <M  ,  H  (I.)> 
ml  — m  — t  — 1  A 


b 

2 


d 

TT 


-jk 


V(y-r  '2‘2 


yc) 


Tm<y>  — rr~r~dy 

(y-yc)  +  d 


Also,  T  (y)  is  approximated  by  four  impulses  as  described  in  Appendix  B. 
m 
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on  the  axis  of  the  slot  (x  *  0,  z  =  0) .  Therefore 


H„(I_)  =  jz  T  (x‘)  (4  +  e  jkRdx’ 

21T  l  n  R3  R2 

wire 


(71) 


Substitution  of  (23)  and  (71)  into  (16)  gives 


T  =  <M  ,  H  (I  )> 
mn  — m  t  — n  A 


_d_ 

2tt 


b 

2 

* 

■  b 


t  (y’> 


wire 


T  (x')  (*\  +  ^)e'jkRdx’dy’ 
n  R3  R2 


Again  Tm(y)  and  T^(x)  are  approximated  by  four  impulses.  The  result  is 

T  =  4  l  l  C^(A )  C*(B)  (^—  +  -&-)  e  JkRmnAB  ^  Ay, 

2  A=1  B*1  m  n  R3  R2  „  2n-2+B  2m-2+A 

mnAB  mnAB 


(m  =  1,2,..  .NA,  n  *=  2.3....NB) 


(72) 


where 


RmnAB  > 


mnAB  ^X2n-2+B  +  (y2m-2+A  “  Yc')  +  d 


(73) 


is  the  distance  from  the  center  of  the  (2n-2+B)th  segment  on  the  wire  to 
the  center  of  the  (2m-2+A)th  segment  on  the  slot.  From  (28)  and  (72),  all 
the  elements  of  [T]  also  can  be  evaluated. 


(D)  Evaluation  of  I 

The  incident  field  considered  here  is  a  plane  wave,  which  is  a 
good  approximation  for  field  radiated  by  a  source  distant  from  the  slot. 
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The  propagation  vector  k  is  assumed  to  lie  in  the  x-z  plane,  as  shown 
in  Fig.  4.  There  are  two  polarizations  for  the  incident  plane  wave. 

One  is  vertical  polarization  for  which  the  electric  field  is  perpen¬ 
dicular  to  the  incident  plane,  or  in  the  y  direction.  The  other  is 
horizontal  polarization  for  which  the  electric  field  is  in  the  inci¬ 
dent  plane,  or  the  magnetic  field  is  in  the  y  direction.  In  our  narrow 
slot  case,  we  may  neglect  the  vertical  polarization  because  it  has  no  y 
component  of  magnetic  field,  which  is  the  only  component  that  couples  to 
region  b.  Hence,  we  only  need  consider  horizontal  polarization. 


Fig.  4.  The  incident  wave  from  region  a. 


For  simplicity,  we  assume  that  a  unit  incident  plane  wave 


propagates  in  the  z-direction,  and  is  of  the  form, 


H1  -  u  e"jkz 

-y 


i  -jkz 

E  -  u  n  e 
—  —x 


(74) 

(75) 


where  n  denotes  the  intrinsic  impedance  of  free  space.  Therefore  the 
tangential  magnetic  field  H1  of  incident  wave  on  the  slot  is 


H 


.  .  =  u  2e 
on  slot  — y 


-jkz 


z=0  ^~y 


(76) 


since  the  slot  is  closed  by  a  conductor,  forming  a  complete  plane. 
Substitution  of  (76)  into  (17)  and  (23)  gives 


I1  =  <M  ,  Hi> 
m  -m  — 


=  2  T  (y)dy 


m 


■  2J,  c>>  4^- 

A=1 


2m-2+A 


where  Cy(A)  and  Ay_  ...  are  defined  as  before, 
m  2m-2+A 


(77) 


IV.  SOLUTION  TO  EQUATIONS  (13)  AND  (20) 

So  far,  we  have  evaluated  all  the  matrices  in  Equations  (13) 
and  (20)  for  our  slot-wire  problem.  The  computer  programs  are  listed 
in  Part  2  of  this  report.  Before  solving  the  slot-wire  problem,  we 
first  reduce  our  problem  to  the  dipole-wire  problem.  In  this  case  we 
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solve  (20)  with v as  known, and  compare  our  moment  solution  with  Kajfez's 
results.  It  is  shown  in  Appendix  C  that  these  two  solutions  agree  well 
with  each  other. 

For  slot-wire  problem,  given  a  set  of  geometric  parameters,  b,  w, 
a,  d  and  y  ,  we  can  solve  equations  (13)  and  (20).  We  must  also  choose 
a  definite  region  (X  <  L/2)  on  the  wire  for  which  the  non-TEM  modes  are 
assumed  to  be  significant.  The  solution  V  determines  the  axial  distri¬ 
bution  of  magnetic  current  on  the  slot  according  to  (8).  The  solution 
I  determines  the  distribution  of  electric  current  on  the  wire  according 
to  (9).  Here  ot^  determines  the  amplitude  of  outgoing  traveling  TEM  cur¬ 
rent  on  the  wire,  and  a^,...  determine  the  evanescent  currents 
on  the  wire.  It  is  found  that  for  the  particular  choice  d/A  <0.1  and 
L/A  >  1.0,  a  change  in  L/A  affects  very  little.  However,  there  is  a 
small  change  in  the  other  a's.  This  indicates  that  the  evanesccr-*'  waves 
are  localized  to  a  region  near  the  slot,  and  their  vat.Mr»s  distant  from 
the  slot  do  not  have  a  measurable  effect  on  the  field. 

V.  EQUIVALENT  CIRCUIT 

It  is  known  that  the  equivalent  circuit  of  a  wire  behind  an  aper¬ 
ture  perforated  conducting  plane  has  the  form  as  shown  in  Fig.  5.  It 
applies  to  the  TEM  mode  at  the  position  x  =  0.  The  equivalent  sources 


Fig.  5.  Equivalent  circuit. 


] 
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and  V^are  due  to  the  incident  wave  from  region  a.  The  two-port 
network  parameters  depend  only  on  the  geometry  of  the  problem.  For 
computational  convenience,  we  will  represent  the  network  in  terms  of 
its  [Y]  parameters. 

In  order  to  determine  the  parameters  of  the  two  port  network, 
we  remove  the  incident  wave  from  region  a,  apply  a  source  to  the  trans¬ 
mission  line  formed  by  the  wire  and  conducting  plane,  and  then  find 
the  relationship  between  the  port  voltages  and  currents.  To  obtain  the 
elements  of  Y  parameters,  we  use  the  definition, 


i  =  T  v  +  Y  v 
1  11  1  *12  2 


i  =  y  V  +  Y  V 
2  21V1  *22V2 


(78) 


(79) 


where  i^,  i^,  and  are  port  currents  and  voltages  at  ports  1  and  2, 
respectively.  The  reference  directions  are  as  shown  in  Fig.  6. 

To  represent  the  excitation,  we  mathematically  apply  a  TEM  wave 
to  the  infinitely  long  transmission  line.  The  incident  fields  are  then 


,T  _  t+t_  -jkx 
-TEM  1  -  e 


E  =  V+e  e'Jkx 
-TEM  -  6 


(80) 

(81) 


where  I+  and  V+  are  the  mode  current  and  voltage,  and  e  and  h  are 


Fig.  6.  Reference  directions  for  i  and  v. 
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transverse  electric  and  magnetic  fields  of  TEM  mode.  These  have  the 
form  [1] 

e  =  2d _ _ 

77  f(y-yc)2  +  (z-d)2]  [(y-yc)2  +  (z+d)2]  ^ 

2  2  2 
d^  +  (y-y  )Z  -  zl 

_  - - - - - - —  „  (a«d)  (82) 

[(y-y  )  +  (z-drj  [(y-y  )  +  (z+d)z] 


“X  ~ 

2d _ z(y-yc} _ u 

v  [ (y-yc) 2  +  (z-d)2]  [(y-yc)2  +  (z+d)2]  ~z 
2  2  2 

d _ d  +  <y~yc>  - 2 _ 

77  [ (y-yc) 2  +  (z-d)2]  [(y-yc)2  +  (z+d)2] 


(83) 


For  the  present  case,  in  which  the  excitation  is  an  incident  TEM  wave, 
every  equation  in  Sec.  II  of  this  report  is  the  same  except  that  (2) 
must  be  changed  to 

H3  =  Ha(M) 

and  (5)  must  be  changed  to 

Hb  =  Hb(-M)  +  Hb(I)  +  H1 

where 

I1  =  =  I+  h  e-jkx  (84) 

->1 

the  final  equations  are  still  (13)  and  (20),  except  for  I  .  Now  I 


becomes 
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I1  =  <-M  ,  H1> 
m  ~m  — t 


b 
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Jb 
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t  (y) 
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~y 


I+  h  dy 
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2 

b 
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T  (y) 
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7T 


2  2 

(y-yc)  +  d 


dy 


I+  r  V  -1  yi“yc  -1  y?yc 

—  I  Cy( A)  [tan  (“V9)-  tan  (^r^)  ]  (85) 

77  A=1  m  d  a 


y  -f  _ 

where  C^(A)  are  defined  as  before,  and  and  y^  denote  the  two  ends  of  the 

ith  segment  Ay^  (i  =  2m-2+A) . 

The  incident  TEM  wave  will  excite  both  outgoing  waves  and  evanescent 
waves.  Because  of  the  assumption  of  a  narrow  slot,  and  the  symmetry  of 
the  problem,  both  outgoing  waves  have  the  same  amplitude. 

We  now  go  through  the  procedures  described  in  Sec.  II  once  again.  The 
only  difference  is  that  the  excitation  vector  now  is  represented  by  (85) 
instead  of  (77).  The  others  are  all  the  same  as  before.  We  obtain  the 
outgoing  wave  by  solving  (13)  and  (20). 

We  turn  to  the  equivalent  network  in  Fig.  6.  The  port  voltage  at 
port  1  is 


v,  =  v  -  a,  Z 
1  1  o 


(86) 


where  Zq  is  the  characteristic  impedance  of  the  transmission  line  formed 
by  the  conducting  plane  and  the  wire,  v*  is  the  mode  voltage  which  we 


applied  to  the  transmission  line  traveling  in  the  +x  direction. 
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i  + 
v  =  I  Z 


The  port  current  at  port  1  is 


At  port  2,  we  have 


vl  =  I  +  ax 


v2  =  v  +  nt1Zo 


i2  =  -  (I  +  a1) 


Because  of  symmetry,  it  follows  that 


Y  =  Y 
11  22 


From  the  reciprocity,  it  follows  that 


Y  =  Y 
12  21 


Hence,  (78)  and  (79)  can  be  reduced  to 


i  =  Y  v  +  Y  v 
X1  11  1  x12  2 


i2  =  Y12V1  +  Y11V2 


Solving  (91)  and  (92)  for  Y^  and  Y12>  we  obtain 


Y11  ’  (11V1  -  I2V2)/(V1  -  V2> 


V12  ’  (12V1  '  1l',2)/(vl  -  V2> 


(87) 

(88) 

(89) 

(90) 


(91) 

(92) 

(93) 

(94) 


Substituting  (86)  through  (90)  into  (93)  and  (94),  we  find 
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where 


+ 

I  +  a. 


11 


2ax  o 


1+  +  a. 


12 


2al  0 


Y  =  1/Z 
o  o 


It  is  clear  from  (95)  and  (96)  that 


(95) 

(96) 


(97) 


Y  =  -  Y  =  Y 
11  12 


where 


+ 

I  +  a. 


Y  =  - 


2a  o 


If  I  equals  unity,  we  have 


Y  "  *  <*  +  ^  T 


(98) 


(99) 


Now,  from  (98),  we  can  plot  the  equivalent  network  shown  in  Fig.  7. 
Finally,  we  evaluate  the  equivalent  sources  in  Fig.  5  when  the  incident 
wave  comes  from  region  a.  Since  the  network  has  the  form  of  Fig.  7, 
we  may  reduce  Fig.  5  to  Fig.  8. 


+ 


«2 


-  > 


-o- 


Fig.  7.  Equivalent  network  at  x  ■  0. 
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Vo  Y 


Fig.  8.  Equivalent  network  for  the  original 
problem. 

i 

If  both  ports  are  terminated  with  matching  loads,  as  shown  in  Fig.  9, 

I 

the  problem  corresponds  to  that  of  an  infinitely  long  wire.  In  this  case,  j 

we  can  solve  matrix  equations  (13)  and  (20)  to  obtain  the  outward  traveling 
current  on  the  wire. 
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From  Fig.  9,  it  is  clear  that 

Vo  -  -  <7  +  2V  <10<» 

where  I  is  the  outgoing  current  when  the  wire  is  infinitely  long  and  the 
wave  comes  from  region  a.  The  Y  of  the  equivalent  circuit  is  given  by 
(99). 


VI.  NUMERICAL  EXAMPLES  AND  DISCUSSION 

A  computer  program  was  written  to  calculate  the  distribution  of 
magnetic  current  on  the  slot  and  electric  current  on  the  wire  for  an 
infinite  wire  case.  The  parameters  of  the  equivalent  circuit  are  also 
calculated.  This  program  is  described  and  listed  in  Part  2.  In  this 
section,  we  give  some  numerical  examples  foT  a  few  cases  of  interest. 

In  Fig.  10,  we  show  the  distribution  of  the  slot  axial  magnetic 
current  along  the  slot  axis  for  the  typical  case  of  b  =  0.5A,  w  ■  0.05A. 

The  incident  plane  wave  is  incident  normally  on  the  slot  from  the 
region  a.  The  infinitely  long  wire  is  located  at  the  position  Yc  m  0 
and  d  *  0.1A.  Also,  the  magnetic  current  for  an  isolated  slot  is  plotted 
in  Fig.  10.  Note  that  the  current  for  the  isolated  slot  is  different  in 
magnitude  from  that  when  the  wire  is  present. 

In  Fig.  11,  we  show  the  distribution  of  the  total  induced  electric 
current  on  the  wire.  The  real  part  is  shown  in  Fig.  11a,  and  the  Imaginary 
part  in  Fig.  lib.  In  each  case  the  total  non-TEM  current  is  also  shown. 
Figure  10  shows  that  the  field  on  the  slot  is  influenced  by  the  presence 
of  the  wire.  Figure  11  shows  that  the  induced  current  on  the  wire  consists 
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lm  (I/XHl) 


Fig.  11(b).  The  imaginary  part  of  current  on  the  wire,  for 
b/A  =  0.5,  w/A  »  0.05,  a/A  ■  0.001,  d/A  -  0.1. 


Fig.  13.  Distribution  of  the  axial  slot  magnetic  current  for 
various  y  /X.  (d/X  =  0.1,  b/X  =  0.5,  w/X  =  0.05, 


Fig.  15.  Equivalent  admittance  for  w/A  =  0.05,  d/A  =  0. 

a/A  *  0.001  (A  =  1  meter,  H  =  1  Ampere/meter) 


Equivalent  source  for  w/A  -  0.05,  d/A  -  0.1, 
a/A  ■  0.001  (A  =  1  meter,  ■  1  Ampere/meter) 
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primarily  of  outgoing  TEM  mode  currents.  The  higher  order  modes  are 
significant  near  the  slot,  but  they  die  out  rapidly,  as  expected.  The 
current  move  than  one-half  wavelength  from  the  slot  is  essentially  a 
single  TEM  mode.  In  Fig.  11,  wL  denotes  the  region  on  the  wire  where  we 
assume  that  the  higher  order  mode  currents  are  considerable.  Changing  wL 
changes  the  total  current  hardly  at  all.  Even  for  the  higher  order  modes 
there  is  very  little  change  as  wL  is  changed. 

In  Fig.  12  and  Fig.  13,  we  show  how  the  position  of  the  wire  with 
respect  to  the  center  of  the  slot  affects  the  distribution  of  magnetic 
current  on  the  slot.  Figure  12  is  for  the  wire  close  to  the  conducting 
plane.  Note  that  a  displacement  of  the  wire  from  the  center  of  the  slot 
changes  both  the  •symmetry  and  the  magnitude  of  magnetic  current.  When 
the  distance  d  is  increased,  it  does  not  change  the  symmetry  much, 
although  the  magnitude  of  magnetic  current  is  still  influenced  by  the  dis¬ 
placement.  This  is  shown  in  Fig.  13. 

The  variation  of  the  equivalent  network  parameters  with  the  dis¬ 
placement  from  the  center  of  the  slot  is  shown  in  Fig.  14  for  d/X  *  0.1, 
b/X  *>  0.5,  w/X  ■  0.05  and  a/X  =  0.001.  As  expected,  they  are  symmetrical 
about  the  centerpoint  of  the  slot. 

Figures  15  and  16  show  the  variation  of  the  equivalent  circuit 
parameters  when  the  length  of  the  slot  is  changed.  For  the  case  w/X  -  0.05, 
d/X  ■  0.1,  a/X  -  0.001,  a  resonance  occurs  in  the  vicinity  of  b/X  **  0.46  when 
the  imaginary  part  of  equivalent  admittance  vanishes,  the  equivalent  source 
reaches  its  maximum.  As  expected,  the  equivalent  admittance  is  inductive 
when  b/X  <  0.46,  and  capacitive  when  b/X  >  0.46. 


Equivalent  source  for  w/X  *  0.05,  b/X  -  0.5, 
a/X  *  0.001.  (X  =  1  meter,  H*  ■  1  Ampere/meter) 


Maximum  power  transferred  to  the  load  for  d/A  *  0.05  and 
d/A  ■  0.10,  and  power  transmitted  through  the  aperture 
when  no  wire  is  present. 
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Figures  17  and  18  show  the  variations  of  equivalent  circuit  parameters 
with  distance  d.  It  is  evident  from  Figs.  17  and  18  that  the  magnitude  of 
equivalent  source  decreases  and  the  magnitude  of  equivalent  admittance  in¬ 
creases  as  the  distance  between  the  plane  and  the  wire  is  increased. 

It  is  interesting  to  compare  the  power  transmitted  from  region  a  to 
region  b  through  the  slot  when  no  wire  is  present  to  the  maximum  power 
coupled  to  the  transmission  line  when  the  wire  is  present.  For  an  isolated 
slot,  the  power  transmitted  from  region  a  to  region  b  is  [6] 

P  =  V[Yns]  V  (101) 

-> 

where  V  is  the  coefficient  vector  of  the  magnetic  current  on  the  isolated 
slot,  and  V  is  its  transpose.  [YhS]  is  the  generalized  admittance  matrix 
of  the  slot  in  half  free  space.  The  asterisks  denote  complex  conjugate. 
Recalling  2[Y^S]V  =  I*,  where  I*  is  the  excitation  vector  given  by  (77), 
we  reduce  (101)  to 

Pt  =  J  V  I1  (102) 

Note  that  P  is  a  complex  power.  The  time-average  power  transmitted  by 

the  slot  is 

Pr  =  Re(Pt)  (103) 

where  Re(Pt)  denotes  the  real  part  of  Pt. 

The  maximum  power  output  in  the  transmission  lines  is  obtained  when 
it  is  terminated  for  maximum  power  transfer,  as  shown  in  Fig.  20.  Here  Z 
is  the  equivalent  impedance 


Z  -  1/Y  *  R  +  jx 


(104) 


Fig.  20.  The  equivalent  circuit  for  maximum 
power  output 


The  loads,  and  Z^*  referred  to  x  =  0,  will  receive  maximum 
power  when 


ZL1  +  ZL2  "  R  “ 
in  which  case  the  power  is 


L  2R 


(105) 


(106) 


the  results  are  plotted  in  Fig.  19.  It  can  be  seen  that  the  power 
transmitted  from  region  a  to  region  b  through  the  slot  can  be  greatly 
increased  by  the  presence  of  a  nearby  wire. 
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Appendix  A 

THE  ELECTROMAGNETIC  FIELD  DUE  TO  TRAVELING  WAVE 
CURRENTS  ON  AN  INFINITELY  LONG  WIRE 

To  evaluate  the  [Z]  and  [T]  matrices,  it  is  necessary  to  know  the 
electric  and  magnetic  fields  due  to  outward  traveling  currents  on  an 
infinitely  long  wire.  The  geometry  and  coordinates  are  shown  in  Fig.  A-l. 
Here  two  traveling  currents  start  from  the  origin  and  travel  to  infinity 
in  the  +z  and  -z  directions.  We  are  interested  in  the  component  of  electric 
field  parallel  to  the  wire,  E^,  and  the  component  of  magnetic  field  perpen¬ 
dicular  to  the  wire,  H^.  So  only  these  two  components  of  the  fields  are 
evaluated.  Because  of  symmetry,  it  is  convenient  to  use  a  cylindrical 
coordinate  system  as  shown  in  Fig.  A-l.  The  point  p(p,<J>,z)  is  the  field 
point  where  the  fields  are  to  be  determined.  For  simplicity,  we  start 


Y 


Fig.  A-l.  Geometry  for  fields  due  to  a  current  along  the 
z  axis. 


R  =  Vp2  +  (z  -  z’)2 


Rl=Vp2  +  (Z  -  L)2 


R  =  Vp2  +  z2 
o 


(A-l) 


(A-2) 


(A-3) 


The  current  is  assumed  to  be  traveling  in  the  +Z  direction  on  the  wire, 


given  by 


I+  =  C  e~jkz  u 


(z  >  0) 


(A-4) 


where  u  is  the  unit  vector  in  the  z  direction. 

— z 

The  expression  for  the  vector  potential  at  point  P  is 


A  =  A  u 

—  Z  2 


u  f  -ikz’e"jkR 
A  =  C  e  J  — - —  dz ' 

z  4tt  r 


-jk(R+z ' ) 


(A-5) 


The  4>  component  of  magnetic  field  at  P  is 

,  9A 

hT  =  ~  (V  x  A)  ,  =  -  — 
9  y  —  0  9p 


If  1  , 

4tt  J  9p  ^  R  )dz 


4 it  J  9p 
0 

L 


..Lf 

4,j  r2  r3)C 


The  integrand  turns  out  to  be  a  perfect  differential.  Integration  there¬ 
fore  gives 
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H+_  C  rpe-J^,)12,=L 

<p  4tt  lR(R  +  z'-z)  , 

z  u 


C£  fe_ 


■jkCR^  +  L)  -jkRc 


4tt 


Vrl  +  l"z)  VVz) 


Cp  r(RL  -  L+Z)e 

A  TT  l  O 


-jko^+z) 


-jkR 

(R  +  z)  e  ° 


4TT  1rl[rl  ~  <L"Z>2]  -  z2) 


o  o 


From  (A-2)  and  (A-3)  we  obtain 


d2  2  d2  2  2 
-  (L-z)  =  Rq  -  z  =  p 


Hence, 


.  r  -jk(R+L)  -jkR 

H*  "  -  4^  [(1  -d  +  f)e  °]  <A~6> 


Now  let  L  -*■  00  and  note  that  lim  —■  =  1.  Now  (A-6)  reduces  to 

L-x»  \ 

+  C  z  “J  kR 

#")e 

O 

This  is  the  <|>  component  of  magnetic  field  at  point  P  due  to  a  current 
traveling  in  the  +Z  direction,  starting  at  z  =  0. 

The  electric  field  can  be  obtained  from  the  magnetic  field  by 
recalling  that 

V  x  H  =  jojeE 

from  which  we  find 

9H 


(A-7) 


Because  of  symmetry  there  is  no  variation  in  the  <\>  direction.  Therefore 

(A-8) 

Substituting  (A-7)  into  (A-8),  we  obtain 


„+  1  1  3  ,  „  . 

EZ  jODE  p  9p  <j)^ 
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z 


C 

4  IT  j  OOP 


o 


+  3> 

R 

o 


-jkR 

o 


(A-9) 


Ez  is  the  tangential  component  of  electric  field  at  point  p  due  to  a 
current  traveling  along  the  +z  direction,  starting  at  z  =  0. 

Now  assume  we  have  an  outward  traveling  current  starting  at  z  =  0 
and  extending  to  z  =  -  °°, 

L  =  C  ejkz  u^  (z  <  0)  (A-10) 

To  obtain  its  fields  and  E^,  all  we  need  to  do  is  to  replace  z  by  -z 
in  (A-7)  and  (A-9).  It  follows  that 


H .  = 


C 

4?rp 


“  -r> 

o 


-jkR 


E 

z 


C  ,jk  jkz  z  . 
(R„  "  R2  '  R3^ 


(A-ll) 

(A-12) 


o  o 

The  total  field  is  the  sum  of  those  due  to  the  +z  traveling  current  and 
the  -z  traveling  current,  or 


where 


E 

z 


+ 

<P  <P 


_C_  -^o 

2-rrp  e 


E+  +  E_ 
z  z 


Ck 

2TTweR 

o 


-jkR 

o 


(A-13) 
(A-14 ) 


P 


and  R  is  given  by  (A-3) . 
o 


(A-15) 
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GENERALIZED  Z-MATRIX  OF  A  WIRE  ABOVE  AN 
INFINITE  CONDUCTING  PLANE 

The  problem  is  to  evaluate  the  elements  of  the  generalized  Z- 
matrix  for  a  finite  thin  wire  above  an  infinite  conducting  plane.  We 
follow  the  procedure  used  by  Chao  and  Strait  [8] .  The  geometry  and 
coordinates  are  shown  in  Fig.  B-l.  For  a  thin  wire,  the  following 
approximations  are  made: 

1.  The  current  is  assumed  to  flow  only  in  the  axial  direction  of  the 
wire. 

2.  The  surface  current  is  assumed  to  be  circumferentially  independent 
and  accounted  for  by  a  total  axial  current  I. 

3.  The  only  boundary  condition  is  that  the  axial  component  of  electric 
field  be  zero  on  the  wire  surface. 


Fig.  B-l.  A  wire  above  a  conducting  plane. 
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With  these  approximations  in  mind,  we  evaluate  the  [Z]  matrix  as  follows. 

Recalling  the  definition  of  Z  ,  we  start  from 

mn 


Z 

mn 


=  -  <1  , 
-m 


W" 


(B-l) 


where  is  the  nth  expansion  function  of  electric  current  on  the  wire, 

and  E  (I  )  is  the  X  component  of  electric  field  due  to  I  .  Here  we  use 
— t  — n  — n 

Galerkin's  method,  so  that  the  testing  function  I  is  the  same  as  expan- 

— m 

sion  function  I  . 

— m 

For  relatively  fast  convergence,  we  choose  triangle  functions  as 
current  expansion  functions,  that  is, 


I  (x)  =  T  (x) 
n 

The  wire  is  divided  into  (2NB  +  2)  segments,  where  NB  is  the  number  of 

triangles  as  well  as  expansion  functions.  Each  triangle  has  unit  height 

with  peak  at  the  point  p„  , . ,  the  direction  of  T  (x)  is  coincident  with 

zn+i  — n 

the  axis  of  the  wire  as  shown  in  Fig.  B-2.  Each  triangle  function  T 


Fig.  B-2.  Triangle  functions  Tn(x),  T^^Cx). 


is  non-zero  only  over  four  segments.  Successive  triangles  overlap  two 

segments  except  at  the  ends  of  the  wire. 

Note  that  ^(1^)  is  evaluated  with  the  conducting  plane  present. 

From  image  theory  [A],  this  is  equivalent  to  the  field  produced  in 

region  z  >  0  by  I  and  its  image  I'  without  the  conducting  plane.  Note 
— n  — n 

that  I'  ■  -  I  as  shown  in  Fig.  B-3.  Therefore  the  x  component  of  elec- 
— n  — n 

trie  field  on  the  surface  of  the  wire  can  be  expressed  by 

3<t> 

E  =  -  1u)  A  -  -r— —  on  the  wire.  (B-2) 

x  nx  dx 

Here  A  ,  d>  are  the  vector  potential  and  the  scalar  potential  due  to 
nx  n 

I  and  its  image  I',  given  by 
— n  — n 


\  »  fr  T  (x’) 

nx  4tt  J  n  v  R 

P„.  , 


rr-)dx* 


L  j-2n+3 

4  ire  J  °n 


2n+3  -jkR  -jkR' 


where 


„  i  d  y ’> 

n  joj  dx' 


R  =  V(x-x')^  +  a^ 


R'=  y/ (x  -  x')2  +  (2d)2 


are  the  distances  from  the  current  I  and  I'  to  the  surface  of  the  wire 

n  n 

respectively. 

Substituting  (B-2)  into  (B-l) ,  we  have 

'■2m+3  , 

Z  -  T  (x)  (ju  A  +  -j-  <j>  )dx  ( 

mn  m  J  nx  dx  Yn  v 


.  d<p  dT 

4-  (T  <J>  )  -  T  ~  +  <P  ~ 

dx  m  n  m  dx  n  dx 


it  follows  that 


2tn+3 

(2m+3 

?2nrt-3 

d(T  d>  ) 
m  n 

*  T 

j  m  n 

+ 

4 

2m-l 

P2m-1 

P 

2m-l 

<J>  dT 
n  m 


The  left  hand  side  of  (B-9)  is  zero,  since  T  is  zero  at  the  ends, 

in 

P,. .  ,  and  ^  •  Hence  (B-8)  can  be  reduced  to 
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mn 


2m+3 


dx 


2n+3 


P  P 

2m-l  2n-l 


n  dT  (x)  dT  (x * )  -jkR  -jkR' 

[joiy  T  (x)T  (x')  +  - ^ - 7— ; — ]  - ~i — =-5— )dx' 

J  m  n  jtoe  dx  dx'  J  4ttR  4ttR'  ' 

(B-10) 


In  evaluating  the  integral  of  (B-10),  triangle  T  is  conveniently  approxi- 

n  dT 

mated  by  four  pulses  as  shown  in  Fig.  B-4,  and  is  also  represented  by 

dT 

four  pulses  as  shown  in  Fig.  B-4.  T^  and  are  approximated  by  four  im¬ 
pulses  as  shown  in  Fi^.  B-5.  The  pulse  amplitudes  are 


Cn(D 


I  Ax2n-1 
Ax2n-1  +  Ax2n 


(B-ll) 


Cn(2) 


Cn(3) 


Ax2n-1  +  1  AX2n 
^n-l  +  ^n 

2  **2n+l  +  Ax2n+2 
Ax2n+1  +  Ax2n+2 


(B-12) 


(B-13) 


Cn(4) 


2  Ax2n+2 


Ax2n+1  +  Ax2n+2 


(B-14) 


D  (1)  =  D  (2)  =  t- - t  ~7~~ 

n  n  Ax2n-1  +  ^X2n 


(B-15) 


D  (3)  -  D  (4)  =  7- 
n  n  Ax 


-1 


2n+l  +  Ax2n+2 


(B-16) 


The  impulse  amplitudes  are  Ax2m-2+i  Cm^i^’  Ax2m-2+i  Dm^  (i=l,2,3,4). 
With  these  approximations,  (B-10)  reduces  to 


Z  “  jcoyZ,  +  - —  Z„ 
mn  J  1  ju)G  2 


(B-17) 


where 
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h  "  J  ^  Cm(A)  Cn(B)  **2«-2*A  ^n-Z+B  ^Zm-Z+A’  Q2n-2+B)  (B-18> 

A*  A  d*x 

Z2  "  E  _^.Dm(A)  VB)  ^m-Z+A  ^Zn-Z+fl  ^Zm-Z+A*  Q2n-2+B)  (B-19) 

A“1  B“1 

Here 

4'(Q1,  Qj)  -  mv  Qj)  "  m±,  Qj)  (B-20) 

and  Green's  function  ^  is  defined  as  [10] 


HQ±,  Qj) 


1 

4,,^ 


(B-21) 


where  Q.^  is  at  the  center  of  the  ith  segment  of  the  wire,  Qj  is  at  the 
center  of  the  jth  segment  of  the  wire,  Qj  is  at  the  point  of  the  image 
of  Qj,  and  R  is  the  distance  from  Q^  to  dx. 


I 
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Appendix  C 

APPLICATION  OF  EQUATIONS  (13)  AND  (20)  TO  THE 
MAGNETIC  DIPOLE  AND  WIRE  PROBLEM 

Figure  C-l  shows  the  problem  to  be  considered  and  defines  the 
coordinates  and  notation  to  be  used.  An  infinite  conducting  plane 
covers  the  entire  z  =  0  plane.  An  infinitely  long  thin  wire  is 
oriented  in  the  direction  of  the  x-axis  at  a  distance  z**d  from  it.  The 
source  is  a  unit  magnetic  dipole  situated  at  the  origin.  What  we  wish 
to  find  is  the  electric  current  on  the  wire.  It  is  easy  to  solve  this 
problem  by  the  method  of  moments,  as  described  in  [5].  The  problem  is 
equivalent  to  a  wire  antenna  scattering  problem  in  which  the  source  is 
an  equivalent  dipole  in  free  space.  The  wire  is  the  scatterer  and 


Fig.  C-l.  Dipole-wire  problem. 
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the  equivalent  dipole  when  imaged  in  the  conducting  plane,  is  2YLl.  Hence, 
we  have  to  solve  matrix  equation 

[Z]  I  =  V  (C-l) 

where  I  is  generalized  current  to  be  determined,  V  is  generalized  voltage 
defined  as 

,  V  *  <W  ,  E*>  (C-2) 

n  n 

where  is  the  nth  testing  function  on  the  wire.  E*  is  the  tangential 
electric  field,  E*,  due  to  2K£.  Finally,  [Z]  is  the  generalized  impedance 
matrix  of  the  wire. 

The  same  analysis  as  used  in  Section  11(A)  applies.  Hence,  we  may 
use  the  same  expansion  and  testing  functions  as  previously,  giving  the  same 
[Z]  matrix  as  expressed  in  (31),  (40)  and  (43).  We  then  solve  (C-l)  to 
obtain  I,  the  distribution  of  the  current  on  the  wire. 

However,  instead  of  using  the  way  described  above,  we  prefer  to 
consider  the  dipole-wire  problem  as  a  special  case  of  slot-wire  problem. 

We  are  then  able  to  check  if  our  theory  is  correct. 

This  can  be  done  by  letting  the  magnetic  current  on  the  slot  be 
the  same  as  the  dipole.  If  the  slot  is  small  enough,  the  magnetic  current 
that  represents  it  looks  like  a  dipole.  Then,  solving  I  from  (20),  we 
obtain  the  solution  of  the  dipole-wire  problem. 

Mathematically,  this  procedure  is  as  follows: 

1.  Let  the  length  of  the  slot  be  very  small,  say  b/X  *  0.001. 
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2.  Let 

M  dy  =  KS.  =  1 

slot 


(C-3) 


where  -M  is  the  axial  magnetic  current  on  the  slot  in  region  b.  Recalling 
(8),  we  may  set 


NA  =  1  (C-4) 

Vx  =  2/b  (C-5) 

for  the  simplest  case.  Equation  (C-5)  results  from  the  fact  that  the 
integral  of  the  triangle  equals  b/2. 

3,  Solve  (20)  for  I,  where  [T ] ,  [Z]  are  defined  as  before. 

Again,  the  coefficient  represents  the  magnitude  of  outgoing  traveling 
current,  and  the  other  a's  represent  the  higher  order  mode  currents  on 
the  wire.  We  compute  the  current  I  from  (20)  for  several  cases.  The  re¬ 
sults  are  given  in  Table  C-l,  along  with  those  obtained  by  using  Kaj fez's 
formulation  [1]. 

Table  C-l.  Magnitude  of  the  outward  traveling  current  induced  on  a 

_2 

wire  by  a  magnetic  dipole.  Units  of  10  /\Kl. 


d /A 

Our  Results 

Kajfez's  Results 

0.05 

-1.1522  +j 0.04 23 

-1.1520 

0.078125 

-0.6727  +j0. 02737 

-0.6721 

0.10 

-0.5011  +j0. 02894 

-0.5007 

0.15 

-0.3085  +j0. 03392 

-0.3100 

0.20 

-0.2209  +j0. 02217 

-0.2214 

0.25 

-0.1664  +J0. 02352 

-0.1707 

Table  C-l  shows  that  our  results  agree  closely  with  Kajfez's  results 
except  for  the  small  imaginary  part.  This  small  imaginary  part  must  be  due 
to  numerical  error  in  our  solution. 
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PART  TWO 


COMPUTER  PROGRAMS 


I  Introduction 


The  program  used  to  compute  the  examples  of  this  report  is  described 
and  listed  in  this  part.  This  program  consists  of  subroutines  ZT,  YH,  MUL, 
SICI,  GAUSS,  TOE,  Subprogram  function  P  and  a  main  program.  Each  is  accompa¬ 
nied  by  an  explanation. 


II  General  Impedance  Matrix 


Subroutine  ZT(NB,  WL,  A,  DD,  Z)  computes  the  generalized  Z-matrix 
for  an  infinite  thin  wire  behind  an  infinite  conducting  plane.  The  input 
variables  are  defined  as 

NB  =  the  number  of  electric  current  expansion  functions  on 
the  wire. 

WL  =  the  length  (in  wavelengths)  where  the  higher  order  mode 
currents  are  considered  to  exist. 

A  *  the  radius  (in  wavelengths)  of  the  wire. 

DD  *  the  distance  (in  wavelengths)  from  the  wire  to  the 
conducting  plane. 


The  output  is  stored  in  two  dimension  array  Z. 
Minimum  allocations  are  given  by 


Complex  Z  (NB,  NB) ,  PP  (N?) 


where 


NP  -  2 (NB  +  1) 


As  described  in  Part  One,  the  electric  current  expansion  functions 
are  divided  into  two  parts.  One  is  the  outward  traveling  currents  which 
travel  from  the  origin  to  infinity.  The  others  are  higher  order  mode 
currents  which  are  assumed  to  exist  in  the  region  |x|  <WL/2.  Therefore  the 
resulting  Z  matrix  is  formed  by  adding  one  column  and  one  row  to  the  Z  matrix 
for  an  open  wire  of  length  WL  behind  a  conducting  plane.  DO  loop  10  com¬ 
putes  this  Z  matrix  and  stores  the  elements  in  Z(I+1,  J+l) .  The  additional 
row  and  column  are  computed  by  DO  loop  20,  and  assigned  as  the  first  row 
and  the  first  column  of  the  resulting  Z  matrix. 

Listing  of  ZT 


SUBROUTINE  ZTCNB.WL.A.DD.Z I 

COMPLEX  ZfSO.SOl .C J.P.PPf 1021 .Zl.Z2.P2 

DIMENSION  CIA). 014) 

COMMON  PI.AK.CJ.C 

L=2*NB 

N*NB-1 

OXsWL/L 

OI  11=0. S/OX 

O(2)*01  1  I 

OI3»=-OU  I 

Of  4 1*0  (31 

AA*2.*00 

K*L-I 

KI*2* ( N8  +  1 I 
OO  3  I  =  1  »K  1 
X*  I-i  . 

3  PPf It*P<K*WL.A. X.O. I-P(K. WL.AA.X.O.) 

OO  10  I*I*N 
00  10  J* 1.1 
Z1*<0..0.) 

Z2*(0. .0.) 

OO  S  K  A*  1.4 
OO  S  KB* 1 .4 
I  A*  2*  I  — 2*KA 
JB*2*J— 24K8 
M*  l+IABS( I A-JBI 
Z1*Z1+C<KA) *C(KB)*PP(M  I 
S  Z2*Z240f KA)*0(KB>*PP(M) 

Z(  !♦!.  J41  l*CJ*DX*OX4I  240*PI*PI*Zl-60.  *Z2> 
Z(34t.I4ll«Z<  !♦  I  •  J4I  » 

CONTINUE 
E«240.*PI*0X 


10 


XN=N+1 .5 
DO  20  Msl.N 
Z1*(0. .0.1 
DO  15  K  A— 1.4 
AM=2*M-2  +  KA 

P2*P(K • ML .A .AM. XN  ) -P(  K  »  ML  •  AA  »  AM  *  XN I 
Z1*Z1*C<KA»*P2 
15  CONTINUE 

ZIM+l  ,  1|*E*21 
20  Zl  l.M+l)aZ<M+l*  1) 

T  =  AKA  A 

CALL  SICXISI.Ct .T) 

Z2=-C1 +CJMSI 
y=ak*aa 

CALL  StCIfSll.C  Il.YI 
Z2=Z2*CI  1-C  J*SI  l 
Z(t. 11-120. *Z 2 
RETURN 
END 


•  I 


III  Generalized  Admittance  Matrix 

The  half  space  generalized  admittance  matrix  of  the  slot  Is  a 
symmetric  Toeplitz  matrix,  due  to  the  symmetry  of  the  slot  and  equal 
segments.  This  means  that  the  Y  matrix  is  of  the  form  as 


[Y]  = 


n-1 


n-1 


Hence,  we  need  only  evaluate  one  column  or  one  row  of  the  [Y]  matrix. 

Subroutine  YH(NA,  B,  W,  Y)  computes  the  first  column  of  the  half¬ 
space  admittance  matrix  for  a  narrow  slot.  The  input  variables  are  defined 
as 

NA  ■  the  number  of  axial  magnetic  current  expansion  functions 
on  the  slot. 


B  ■  the  length  (in  wavelengths)  of  the  slot. 
W  ■  the  width  (in  wavelengths)  of  the  slot. 
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The  output  is  stored  in  one  dimension  array  Y. 

Minimum  allocations  are  given  by 
Complex  Y(NA),  PP(NY) 

where 

NY  =  2(NA  +  1) 

DO  loop  3  is  used  to  compute  the  t p  functions  and  stores  them  in  PP  so  that 
repeated  computation  can  be  avoided. 

Listing  of  YH. 


SUBROUTINE  YHf  NA.R.Y.Y) 

COMPLEX  YC50» .C J.P.PPI 1021 .Yt .Y2 
DIMENSION  C(4).D(4> 

COMMON  PI.AK.C J  .C 
A=k/4. 

K=2*NA+1 

DV=B/IK+1) 

D<1 1=0  .5/CY 
0«2)=0(  1) 

D«3I=-DI  1  > 

0(4t=D(3> 

KI=2*{ NA41 I 
DO  3  I*I.K1 
X=  1  •*<  l-l» 

3  PP< I»=PC K.B.A.X.O. I 

DO  10  1=1. NA 

Y I = (O • .0. ) 

Y2=<0. .0.  ) 

DO  5  K  A=  1.4 
OO  S  KB=1 .4 
I  A=  2*  I  — 2*K  A 
M=l«-IAfiS(  IA-KBI 
Y1«YUC(KAI  *CIKB»*PP<MI 

5  Y2=V240 I K  A| *0 (K  E I *PP<  M  t 

10  YC  I  |=CJ*DV*DV*( VI/30.-V2/C  120.*PI*PI» ) 

RETURN 
ENO 


63 


IV.  Subroutine  TOE 


For  a  symmetric  Toeplitz  matrix  we  have  a  recursive  system  of 
equations  [10]  as  follows: 


=  Y-1  Y 
*1  1  2 


A(1)  =  Y'1 

A(m+1)  =  [1  -  C0(m>>2 j_1A(m) 
m 


m 


♦isl)  ■  -A<m+1)  <  J  ♦«  Vs«  -  w 

J  X 


,  (m+1)  (m)  (m)  (m+1)  ,  . 

<p\.  '  =  -  <]'  .  i  ip  . '  1  <  r  <  m 

yr+l  r+1  ra-r+1  m+1  —  — 


Then 

the  NxN  matrix 

[Z]  =  [ 

Z11 

Z12  *  ’ 

’  Z1N 

Z  2 1 

Z  2  2  ‘  * 

Z2N 

ZN1 

7N2  -  ‘ 

’  ZNN 

where  the  elements  of  [Z]  are  given  by 
i  (N) 


Z11  =  A 


'rl 


-  •  a<« 

r-1 


2  <  r  <  N-l 


z  =  z  .  .  +  cr'  a 

rs  r-l,s-l  r-1 


(N-l)  A (N)  ,„<N-])_  lh(N-l)  A(»)  ^(N-l) 


'Vl  ^N-r+1 


KN-S+1 


2  <  r. 


S  <  N 


64 


Note  that  the  inverse  matrix  [Z]  of  [Y]  is  also  symmetrical,  so  we  need 
only  compute  the  elements  of  the  lower  left  triangle,  and  then  set  the 
upper  right  triangle  equal  to  the  transpose  of  the  lower  left  triangle. 

Subroutine  TOE  (N,  Y,  Z)  is  used  to  invert  a  symmetrical  Toeplitz 
matrix  [Y] .  The  input  variables  are  defined  as 


N  =  the  order  of  the  matrix  to  be  inverted 
Y(N)  =  the  first  column  of  the  matrix  to  be  inverted. 

The  output  is  the  inverse  matrix  of  [Y]  and  is  stored  in  two  dimension 
array  Z(N,N) . 

Minimum  allocations  are  given  by 
Complex  Z(N,N) ,  Y (N) ,  PS(N,N),  DL(N) 

DO  loop  40  evaluates  the  half  elements  of  [Z],  except  for  Z^,  which  is 
evaluated  by  statement  Z(l,l)  =  DL(N) .  The  other  half  elements  are  given 
by  the  statement  Z(J,I)  =  Z(I,J). 

Listing  of  TOE. 


SUBROUTINE  TQEIN.Y.ZI 

c  THIS  IS  YHE  PROGRAM  TO  INVERT  A  TOEPLITZ  MATRIX  YINXNl  INTO  ZCNXNI 

COMPLEX  Z  (50*50  )  •  Y  (  50  »  .  A,  PS<  50 . 50  >  . 0L  C  50 1 
Nl=N-l 
N2=N1- 1 
DL«  I»=l  ./Y<  I) 

PSII.l  »*CL( 1 ) *Y  <21 
00  20  M* 1  »N2 
Ml*M+l 

OLIM1  »  *0LIM>/U  .-PS<M,M)*PSCM.M») 

A*<0..0.  > 

OO  10  K«  1  »  M 

|0  A=A*PSIM.K >*Y(M-K42I 

PSIMI  .Ml  I=-0L<M1  »*<  A-VIM4  2M 
00  30  K*t.M 

30  P SI  MI . K»*PSCM,K»-PSIM.M1-K»*PSIM1  .Mil 

20  CONTINUE 

OL  (  Nl  =  OL  <  M  I/O  .-PS  INI  ,Nl  )*PS(N1  ,Nl  )» 

Z<  1. 1 >  =DL INI 
00  40  I *2  »N 
11*1-1 
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Z( I  ,  I >=-PS(Nl , I l ) *  DC (N I 
ZC  l.I»*Z(t.l) 

00  40  J=2. I 

Z<  I  .  J>=Z(  1 1  ,J-1  )*PS<Ni  ,  u  )*DL(N)*PS(N1 ,  J-t  |  -PS(N  1.N-I4  1  )*0L(NI 
1  4PSIN1 .N-J+l ) 

Z< J,I)sZ«I» Jl 

RETURN 

ENO 


V.  Subroutine  GAUSS 


Subroutine  GAUSS  (N,  A,  B,  EPS,  ISW)  is  used  to  solve  a  linear 
system  of  equations  with  complex  coefficients  by  the  method  of  Gaussian 
elimination.  The  coefficients  of  unknowns  are  stored  in  two  dimension 
array  A,  and  the  constants  are  stored  in  one  dimension  array  B.  Again 
the  solution  of  the  equations  are  stored  in  B.  N  is  the  number  of 
unknowns  as  well  as  the  equations.  EPS  is  a  small  constant  and  ISW  is  a 
message.  If  the  absolute  value  of  the  pivot  of  column  is  larger  than 


EPS,  ISW  equals  to  1.  Otherwise  ISW  =  0. 

Minimum  allocations  are  given  by 
Complex  A(N, N)  B(N) 

Listing  of  GAUSS 

SUBROUTINE  GAUSS! N, A, B *  EPS t ISW » 

C  TO  SOLVE  EQUATIONS  AX=8 

COMPLEX  4(50.50  UR(SO)  .C.T 
NM  1=N-  1 
OO  10  K= 1 »NM1 
C*(  0»  *0* } 

00  2  I=K.N 

IF (CABS < A( I «K)  ) »LE • CAB SIC ) I  GO  TO  2 
C=  A  <  l  .  K  » 

10*1 

2  CONTINUE 
1F(CABS(CUGE.EPSI  GO  TO  3 
IS  4*0 

RETURN 

3  IF(IOaEO.K)  Gb  TO  6 

00  4  J=K  »N 

T*A(Ki J) 

A(K.J)*A(  10. J) 

4  A((0.J)*T 


1 


1 
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T*B(K I 
0<  K) =8  < 10) 

0( 1 0) =  T 
6  KP  1  =K  ♦  1 

C=l./C 
0<K)=B<K)*C 
00  10  JsKPl .N 
A(K.J)=A<K. J)*C 
00  20  I  =  KP1 ,N 

20  A<  I.J)=A( I. J)-A< I.K)«A<K. J) 

10  01 J)*B( Jl-AIJ.K  )*B(K) 

81 NI=B(N)/A(N ,N> 

00  40  Ksl.NMl 
I  =N-K 
C=<0..0.) 

I  p  I =  I  ♦  1 
00  SO  J—  IP  1  »N 
50  C=C4A( I . J)*e( Jl 

40  BC I )  =  8  < 1  1— C 

IS*=1 
RETURN 
END 


VI  Subroutine  MUL 


Subroutine  MUL(L,  M,  N,  A,  B,  C)  is  used  to  multiply  matrix  A 
by  matrix  B  and  stores  the  result  in  C.  Here  A  is  an  (L  *  M)  array, 

B  is  an  (M x N)  array  and  C  is  an  (L * N)  array. 

Minimum  allocations  are  given  by 
Complex  A(l,,  M)  ,  B(M,  N) ,  C(L,  N) 

Listing  of  MUL. 

SUBROUTINE  MUL1L.M.N. A.B.C  ♦ 

C  TO  00  MULTIPLICATION  C«A*B 

C0MR.E  X  A  (50*  50  )»BI50»50I»C<50  »50  I  •• 

00  20  I  *  1  »L 
OO  20  K*1*N 
H*IO*iO») 

00  10  Jsl«M 
»*A(I # JI4BI J.KI4W 
C( 

RETURN 
ENO 
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VII  Subrotulne  CISI 


Subroutine  CISI  (SI,  Cl,  X)  is  used  to  compute  the  sine  and  cosine 
integrals. 

00 

* 

^  v  COS  V  J 

Ci(x)=-  — —  dv 

-X 

X 

Si  (x)  =  I  — - —  dv 

0 


This  subroutine  is  described  in  [11]. 

Listing  of  CIST. 

SUBROUTINE  SICI CSI.CI.XI 
Z=ABS( XI 
IEIZ-4.1  1.1.4 
Vs<4.-ZI*<4.«-ZI 

SI=X*I ( I (( 1 .753  141E-94V41 . 568988E- 7 >* Y4 1 . 3741 68E-5I 4V46.93 9889 

1  E-4l*Y*1.964882E-2 l*Y 44. 395509E- 1  I 

CI»It  5.7721  56E-144LCG(ZJI/Z-ZM  CC  (  1 1 . 386985E-  1  0* V4  1  .584  996E-8) 4V 

2  ♦1.725752E-6I4Y41 .  185  999E-4  )  *Y*4 .99092E-3I  *Y*  l  .  31  S308E-1 1  »*Z 
RETURN 

SI^SINIZ) 

V«COS( Z I 
Z*4./Z 

U*M((C<<<4 .048069E— 34Z-2 .2791 43E—2 l*Z+5 .51507E— 21 4Z—7 • 26 1 646E— 2  ) 
O  *Z+4.987716E-2)*Z-3.332519E-3I*Z-2.314617E-2»*Z-1.134958E-5»4Z 

4  4  6.25001  IE— 2I4Z42  •  5839  89 E -10 

¥=(<((((((< —5.1 08699E- 34Z4 2. 8191 79E-2 I ♦ Z-6. 537283E— 21 *Z 

5  47.902034E-2 |*Z-4 .4004 16E-2 I 4Z-7 . 945556E-3 I4Z42 .60 1293E-2  )«Z 
L  -3.764E-4I *Z-3 • 1 2241 8E-2 I *Z- 6. 646441 E-7) *Z42 .5E-1 

CI*Z*1 SI4V-Y4UI 

SI»-Z*f  SI*U4V*V)41 .570796 

RETURN 

END 


An  accurate  evaluation  of  the  scalar  Green's  function  of  (B-21) 


is  developed  in  [9]. 

Subprogram  complex  function  P(N,  B,  A,  X,  Y)  computes  the 
function 


ip  = 


1 

4ttA£ 


X+AS./2 

X-A£/2 


dx' 


12  2 

where  A£  =  b/(nfl)  is  the  integral  length  and  R  =  ya  +  (x'-y)  is  the 
distance  from  source  point  dx'  to  the  field  point.  This  is  shown  in  the 
following  figure. 


With  reference  to  the  figure,  the  input  variables  are  defined  as. 

N  *  n,  an  integer 

B  ■  b,  a  length  (in  wavelengths)  which  gives  the  subsection 
interval  A£  =  b/(N+l). 


A  ■  a,  the  transverse  coordinate  of  the  field  point. 
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X 

Y 


x, 

y» 


the  longitudinal  coordinate  of 
the  longitudinal  coordinate  of 


the  source  point, 
the  field  point. 


Note  that  the  integrand  is  expanded  in  a  Taylor  series  for  R  <  5A£,  and 
in  a  Maclaurin  series  for  R  _>  5A£.  This  is  done  by  the  statement 


IF  (RR.  GE.  10)  GO  TO  100. 


Listing  of  subprogram  function  P. 

COMPLEX  FUNCTION  P(N,B.A,X,Y> 

C  THIS  IS  THE  PROGRAM  TO  CALCULATE  THE  INTEGRAL  EXPC-JKRI/R 

CONPLEX  T.J.TT  .CEXP 
PI=3.141S92C5 
AK  =  2.*PI 

D=E/( N ♦ 1  I 
AL=0/2. 

DZ=D*ABS( X-Y> 

RsSQRTI A*A+OZ*D Z> 

RRsR/AL 

IF  <  RR*  GE  *10  1  GO  TO  100 
G=SORT(A*A*<OZ-AL>*<OZ-AL)  » 

H  =  SQRT  (  A*  A*  (DZ+  ALt  *  f  OZ  +  AL  •  > 

IFIOZ.GE.ALI  GO  TO  9 

BI  =  ALOGI  (  <DZ4AL*HI*CG*-AL-DZ»  >/t  A* A  »  > 

GO  TO  99 

9  Bl=ALOGC 10Z*AL4H)/C0Z-AL*G»» 

99  B3=CDZ4AL l/2.*H'H  AL-DZ  >/2.*G*A* A*BI/2 • 

B 4=2  .  *  AL  * A* A*  <  2  *AL  **346 .*AL*D  Z  *DZ  >  /3. 

T=CEXP C  —  J*AK*RI /C8.*Pl*AL  » 

V=AK**3/6.*(B4- 3.*R*B3+3. *R*R *D-R * *3*B 1  I 

P=T*CBI-J*AK*  CD-R* B I l-AK* AK/2 •* 183-2. *R*0*R*R*BI )*J*V» 

GO  TO  200 

100  C= AL/R 
C  2=C*C 

e»oz/r 

E2=E*E 

E4=E2*E2 

A l*C/6 • * 1 3. *E 2— 1>4C*C2*C3. -30 .*E2*35. *E4 ) /40 • 

A0»  I  •  *C*A 1 

A2=-E2 /6  •— C  2/40  •  *1  t-1  2.  *E2+t  5.  *E4  > 

A3=C/60.*(3.*E2-S.*E4 I 
A4=E4/ 120. 

BK=AK* AL 
BK2=BK  *BK 

TT=CE  XP  C- J*  AK*R  I  /<  4  .*P  I  *R  I 

P=TT*( A0*  J*BK*A 1+BK24A2+J *BK2*8K*A 3*BK 2*BK2*A 4  I 
200  RETURN 
ENC 
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IX  Main  Program 

The  main  program  computes  the  complex  coefficients,  which  determine 
the  axial  magnetic  current  on  the  slot,  which  determine  the  electric  cur¬ 
rents  on  the  wire,  and  the  equivalent  circuit  parameters  of  the  transmission 
line.  The  main  program  calls  the  subroutines  ZT,  YH,  TOE,  MUL,  GAUSS  and  CISI, 
which  are  described  previously  in  Part  2. 

One  data  card  is  read  in  the  main  program  according  to 
READ  (1,8)  NA,  NB,  B,  WL,  W 
8  FORMAT  (2X2,  3F4.3) 

The  input  data  are  defined  as 

NA  =  the  number  of  magnetic  current  expansion  functions. 

NB  =  the  number  of  electric  current  expansion  functions. 

B  *  the  length  of  the  slot  (in  wavelengths). 

WL  =  the  length  (in  wavelengths)  near  the  slot  on  the  wire 

where  the  higher  order  mode  currents  are  considered  to 
exist . 

W  ■  the  width  of  the  slot  (in  wavelengths). 

Another  input,  the  distance  D,  is  assigned  by  DO  loop  200. 

Minimum  allocations  are  given  by 
Complex  T(NA,  NB) ,  TT(NB,  NA) ,  Z(NB,  NB) , 

YZ(NA,NA),  GZ(NB,  NB) ,  YT(NA,  NB) , 

VM(NA),  VI (NB),  Y (NA) ,  YTI(NA) 


DO  loop  40  computes  the  [T]  and  [f]  matrices.  DO  loop  60 
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Is  used  to  reduce  the  order  of  matrix  equations,  NB  by  NB  to  (— j — )  by 

(— |— ) ,  if  NB  is  odd,  or  (—  +  1)  by  +1),  if  NB  is  even.  This  is 
due  to  the  symmetry  of  the  problem.  It  is  done  to  save  computation. 

DO  loop  120  computes  the  distribution  of  currents  twice.  The  first  time 
is  for  incidence  from  region  b,  the  incident  wave  being  a  TEM  traveling 
wave.  The  second  time  is  for  incidence  from  region  a,  the  incident  wave 
being  a  plane  wave.  The  current  solutions  are  printed  out  each  time 
because  they  are  not  kept  in  permanent  storage. 

In  this  main  program,  we  compute  only  the  case  ycM  =  0.  This 
is  done  by  assign  statement  early 

YC  =>  0 


Listing  of  main  program. 

•  JOB  NA IHENG 

C  THE  MAIN  PROGRAM  TO  SOLVE  SLOT-WIRE  PROBLEM 

COMPLE  X  T!50. 50  I • TTI50 .50  ) ,ZI 50 ,50  » ,V  Z! 50 ,50) . GZ«50 .50  I 

1  .  V  1(501.  VS  (50)  »  VM  (50).VI(50).Y(  50  )  «C  J  .  E  I  ,  E2.  H.CEXP  *YE  . 

2  OGZC  50.501  .YTI  (50)  .GG<50  .SO),yT!SO.SO  ).  V0.AM2) 

DIMENSION  CIA  ) 

EQUIVALENCE  (OGZ.TI  .1  GG.GZ) 

COMMON  PI.AK.CJ.C 
PI-3.14159265 
AK-24P1 
C J-CO. ,t.) 

A*>  0.00  l 
VC-0. 

READ! 1,8)  NA.NB.B.WL.W 

8  FORMAT  (2I2.3F4.3) 

DO  200  NN— 1  .2 
0—0 .1  ANN 

WRITE! 3.9)  NA.B.W 

9  FORMAT  I  //IX  .•NA— •  .  I  2.  3X  ,•  B«»  .  F6.4 , 3X.  •  W»*  .F6.4) 

WRITE! 3. 10)  NB.WL.A.D 

10  FORMAT!  I X.«NB*«  .1  2  •  3X  .  *  XL—*  .F6.4.3X.'  A—*  ,  F6.4 , 3X  .*  D— *  •  F6.4  I 

C  A-THE  RAOIUS  OF  THE  WIRE 

C  W«T!C  WIDTH  OF  THE  SLOT 

C  NA-THE  NUMBER  OF  EXPANSION  FUNCTIONS  ON  THE  SLOT 

C  NB-THE  NUMBER  OF  EXPANSION  FUNCTIONS  ON  THE  WIRE 

C  B-THE  LENGTH  OF  THE  SLOT 

C  0— THE  01STANCE  BETWEEN  THE  SCREEN  AND  THE  WIRE 
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C( I >=0.25 
C( 21=0. 75 
C(3 »=C(2» 

C  C  4  >=C  (  1  I 
0X=ML/(24NB> 

OY*B/( 24NA42I 

CALL  ZT( NB . fcL.A.D.Z) 

CALL  YH(NA.B.W.V) 

C  Y=A0MI TT ANCE  MATRIX  OF  A  SLOT  IN  HALF  SPACE 

DO  12  t  =  1  . N A 
12  VI  I  1=2 •  4V 1 1  I 

CALL  TOEINA.V.VZ) 

C  TO  COMPOTE  COUPLE  MATRICES.T  AND  TT 

IS  WRITE! 3. 20)  YC 

20  FORMAT I///3  X. * YC-*  *F4  »  2  > 

OO  40  M=1,NA 
N1=NB- 1 
OO  30  Nsl.Nl 
E1=(0. .0. > 

E2=C0. .0.) 

C  VS=THE  EXCITATION  VECTOR 

VSI M)=(0..0 .1 
OO  2S  K  A=  1.4 
MA=2*M-2«KA 

VV=(24MA—  1 1  *OY/?,«H  -  -VC 

YP*<YY40Y/2  .1/0 

YN= ( YY— 0V/2 . ) /O 

AT* AT  ANIYP)— ATANI YN I 

VS(M)=VS(M>  —CIR  A  >  4  AT/P  I 

RY*SQRTI  YV4VY+D4D> 

H*CEXPI  — CJ*  AK4RVI 
E2*E2*C(KA)  4M/I  RV4RV1 
OO  25  KB*  1.4 
NKB*24N-24KB 
X*!24NKB-1 140X/2.-WL/2. 

R*SQRT I X4X4VV4VV4DA01 

E1=E1*C(KA» 4C(KB>4(  l/R 4434CJ 4AK/I  P  *R»  )  4CEXPI-CJ4  AK4R» 
25  CONTINUE 

TIM.N4I »=D4DX4DV4EI/AK 
TTIN41 ,M»*-TCM.N41 I 
30  CONTINUE 

TIM.l ) *040Y4E2/PI 
TTI l.M )*— Tl M»  1) 

40  CONTINUE 

CALL  MULINA.NA.NB. VZ.T. YT » 

CALL  MUL (NB  »NA. N8.TT • YT  *GZ  ) 

I N*NB/2 

L*IN4l 

OO  60  t*l.L 

OGZI I.  1 »*GZ(I.l l-ZC  I. 1 1 
00  60  J«2.L 
J1*NB- J42 


DGZII  .  J)=GZII  .Jl-Z!  It  J  I  *GZI  I  •  J1  l-ZII.Jl  1 
IF! 2*1 N. EG. NO )  GO  TC  SS 
GO  TO  60 

55  OGZI  I  •  L  )  =GZ  d.L)— Z((.L) 

60  CONTINUE 

00  120  K*  1 1 2 
00  65  1=  1,L 
DO  65  J= 1  .L 
GG ( I.J  >=DGZ! I.JI 
65  CONTINUE 

IF(K.EQ.l)  GO  TO  85 
WRITE ( 3  t 130 1 
00  75  1  =  1  .NA 
75  VSII)=4.#DY 

GO  TO  90 

85  WRITE! 3# 125  I 

90  CONTINUE 

OO  95  1=1 .NA 
YI(I)=(OttO.) 

DO  95  J=  1  <NA 

95  YIII)=YZII.J)*VSIJ»4YIIII 

OO  100  1=1#  Nfl 

vi  m=  (o#  .0.) 

DO  100  J=1.NA 

100  VI  m=TTfl,  J»*YIC  J)«I(  II 

CALL  GAUSS! L. GG  .VI » IE— 11. 1 SW 1 
IF!  ISW .EG . 1  )  GO  TO  102 

MR I TE (3.101  I 

101  FORMAT !  •  ISW=0«> 

STCP 

102  WR I TE ( 3 . 1 35 1 

00  103  1=2. L 

103  VI! N0— I42)=VIII ) 

WRI TE ( 3 .  140  I  VIII) 

WRITEI3. 1451 

WP  I  TE  (  3  •  1 60  M  V I  <  I  1  .  1=  2  .  N0  I 
OO  105  1=1. NA 
vTim  =  <o..o.) 

00  105  J  = 1 • N8 

105  YTUI  »=YT(  I  .J)*VI(  Jl+YTH  X  » 

OO  110  1=1. NA 
VMI  I)=VI(I)-YTI  II) 

VM!  I  >  =  VMI  1 1/376.99 

110  CONTINUE 

WRITE! 3. ISO  I 

WRITEI  3.  160)!  VMI  I)  .1=1.  NA  > 

120  AI IK) =VI II) 

Z0=60. 4AL0GI 2 .*0/A ) 

VE*-I1 41/4111 ) ) /I2 «*Z0 ) 

VO=-AI ! 2)*I Z0*1 /! 2 • *VE» > 

V0=2*V0 
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WRITE! 3. 165) 

WRITE (3*1 70  )  YE. VO 
REYsREAL! l./YEI 
RV=CABSIVO> 

PL=RVWRV/!2*REY I 
WRITE! 3.  180  I  PL 
200  CONTINUE 

1 25  FORMAT!///*  INCIDENT  WAVE  IS  TEM  WAVE*) 

130  FORMAT (///•  INCIDENT  WAVE  IS  PLANE  WAVE*) 

135  FORMAT!/*  ELECTRIC  CURRENT  ON  THE  WIRE*) 

140  FORMAT 1/3X. *THE  OUTGOING  TRAVELING  CURRENT® */3X. 2F 16.6) 
145  FORMAT !  /  3X. *NON-TE M  CURRENTS  ARE*) 

150  FORMAT  !/•  AXIAL  MAGNETIC  CURRENT  ON  THE  SLOT®*) 

160  FORMAT! 3X.2F16.6) 

165  FORMAT!//*  THE  EQUIVALENT  NETWORK  PARAMETERS*) 

170  FORMAT  !  4X.  *V®*,2El  6. 6/3  X,  •  V0  =  *  .2016.61 
160  FORMAT!  •  THE  POWER  TRANSMITTED  P=*.E!5.6> 

STOP 

END 


(DATA 


N A® 2 1  6  =  0  . 5000  W=0.0500 

NB®32  WL=1  •  6000  A=0  .00  10  0=0.1000 


VC*0 .00 


INCIDENT  WAVE  IS  TEM  WAVE 
ELECTRIC  CURRENT  ON  THE  WIRE 

THE  OUTGOING  TRAVELING  CURRENT® 


-0.20536* 

0.080070 

NON-TEM  CURRENTS  ARE 

-0.000876 

-0.004263 

—0.00  1889 

-0.007896 

—0.002893 

-0.01 1146 

-0.003648 

-0.013871 

-0.003880 

-0.015925 

-0.003312 

-0.017164 

-0.001687 

-0.017454 

0.001205 

-0.016653 

0.00S5I 1 

-0.014597 

0.0  11303 

-0.01 1058 

0.0  16590 

-0.005676 

0.027366 

0.002187 

0.037724 

0.01 3782 

I 
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0.0=01 10 
0.066235 
0.062420 
0.066235 
0.0501 19 
0.037724 
0.027366 

o.o  lesoo 
0.011303 
0 .00551 1 
O.OC  1205 
-0.001687 
-0.003312 
-0.003880 
-0.003648 
-0.002893 
-0.001889 
-0.000876 


0.031745 
0.062589 
0.  100017 
0.062589 
0.03  1745 
0.01 3782 
0.002187 
-0.005676 
-0.01 1058 
- 0 . 0 l 4597 
-0.016653 
-0.017454 
-0.01 7164 
-0.015925 
-0.013871 
-0.01  1146 
-0.007896 
-0.004263 


AXIAL 


MAGNETIC  CURRENT  CN  THE  SLOT 


-0.166122 

-0.224527 

-0.299857 

-0.360378 

-0.415273 

-0.461746 

-0.500140 

-0.529987 

-0.551210 

—  0 . 5  6  383  8 
-0.568021 

—  O  .56  383  8 
-0.551210 
-0.525987 
-0.500140 
-0.461746 
-0.4  15272 
-0. 360378 
-0.255856 
-0.224527 
-0.  166121 


0.071415 
0.092303 
0.  118933 
0.  137290 
0.  151804 
0.  161629 
0. 167483 
0.  169945 
0. 170065 
0.  169204 
0.  168718 
0. 169204 
0. 170065 
0.  169945 
0. 167483 
0. 16 1629 
0. 15  1804 
0. 137289 
0. 118933 
0.092303 
0.071415 


i 


INCIDENT  VAVE  IS  PLANE  4 AVE 

ELECTRIC  CURRENT  CN  THE  HIRE 

THE  OUTGOING  TRAVELING  CURRENT* 
0.225198  -0.091860 
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NON-TEM  CURRENTS  APE 
0.001035 
0*002198 
0*003329 
0.004162 
0*004397 
0*003732 
0.001890 
-0*001353 
-0.006148 
-0.0  12563 
-0.020602 
-0.020252 
-0.041627 
-0.055279 
-0.073243 
—  0.05  1586 
-0.072243 
-0.055279 
-0.041627 
-0.020252 
-0.020602 
-0.012563 
-0.006148 
-0.001353 
0.001890 

0.002732 

0.004397 

0.004162 

0.003329 

0.002198 

0.001035 


0.00471  4 
0.008733 
0.01 2325 
0.015331 
0.01 7587 
0.018932 
0.019216 
0.018285 
0.015959 
0.01 1992 
0.005995 
-0.002728 
-0.015538 
-0.03531  8 
-0.  069199 
-0.110264 
-0.069199 
-0.03531 8 
-0.015538 
-0.00 2728 
0.005995 
0.01  1992 
0.01 5959 
0.018285 
0.01 9216 

0.018932 
0.01 7587 
0.015331 
0.012325 
0.008733 
0.004714 


AXIAL  MAGNETIC  CURRENT  ON  THE  SLOT= 
0.182180  -0.054367 

0.246201  -0.076092 

0.328772  -0.  103953 

0.395092  -0.127723 

0.455233  -0.150033 

0.506132  -0.169758 

0.548168  -0.186742 

0. 560834  -0.200546 

0.604052  -0.210811 

0.617862  -0.217170 

0.622435  -0.219330 

0.617862  -0.217170 

0.604052  -0.210812 

0 • 5808  34  -0.200546 

0.546168  -0.186742 


0.506131  -0.169758 
0.455232  -0.150033 
0.395091  -0.127723 
0.328771  -0.103953 
0.246200  -0.076092 
0.ie?179  -0.054367 


THE  EQUIVALENT  NETWORK  PARAMETERS 

V=  0.507528E-02  0.259205E-02 

V0=  -0.171041E  03  0.907327E  02 

THE  POWER  TRAN  EMIT  TED  P  =  0.1I9943E  03 


NA=2 1  8=0.5000  W=0.0500 

N8=3 2  WL= i * (000  A  =  0  .00  10  0=0.2000 


YC=0.00 


INCIOENT  WAVE  IS  TEM  WAVE 
ELECTRIC  CURRENT  CN  THE  WIRE 


THE  OUTGOING  TRAVELING  CURRENT* 


-0.08  11 29 

0.060615 

NON— TEM  CURRENTS  ARE 

-0.00127* 

-0.00 1675 

-0.002684 

-0.003272 

-0.004092 

-0.004987 

—0 .006232 

-0.006844 

— 0.005788 

-0.00  8815 

—0.005431 

-0.010814 

-0.003840 

-0.01 2685 

-0.000731 

-0.014189 

0.004119 

-0.014994 

0.0  10866 

-0.01 4665 

0.0  19595 

-0.01 2641 

0.030348 

-0.00821  1 

0.0421 54 

-0.00  0484 

0.05e049 

0.01  1584 

0.0 15240 

0.029302 

0.089693 

0.046771 

0.0  75240 

0.029302 

0.058049 

0.01 1584 

0.043154 

-0.000484 

0.030348 

-0.00821  1 

0.019595 

-0.01 2641 

0.010866 

-0.014665 

0.004119 

-0.014994 

-0.000731 

-0.014189 

-0.003840 

-0.01 2685 

-0.005431 

-0.010814 

-0.005788 

-0.008815 

-0  .005232 

-0.006844 

-0.004092 

-0.004987 

-0.002684 

-0.003272 

-0.001276 

-0.00 1675 

AXIAL  MAGNETIC  CURRENT 

ON  THE  SLOT 

-0. I 13120 

0.069438 

-0. 153494 

0.092154 

-0.205619 

0.121465 

-0. 247950 

0.  144032 

-0.286674 

0.  163980 

-0.319843 

0.  180350 

-0.347610 

0.  193531 

-0.369537 

0  .  20  355  5 

-0.385399 

0. 210579 

-0. 394995 

0.214727 

-0.398207 

0. 216098 

-0.394995 

0.214727 

-0. 385399 

0.210579 

-0,369536 

0.  20  3555 

-0.347609 

0.  193531 

-0.3  19842 

0. 180350 

-0.286673 

0. 163980 

-0.247949 

0.  144032 

-0. 20561 9 

0. 121465 

-0.153493 

0.092154 

-0. 1  13120 

0.069438 

INCIDENT 

WAVE  IS  PLANE 

WAVE 

ELECTR  IC 

CURRENT  ON  THE 

WIRE 

THE  OUTGOING  TRAVELING  CURRENTS 

0.1 28560 

-0.096623 

NON-TEM  CURRENTS  ARE 

0.002036 

0.00  2664 

0.004275 

0.005203 

0.006SI 0 

0.007930 

0.008312 

0.010379 

0.009185 

0.014005 

0. 008605 

0.01 7170 

0.006064 

0.020127 

0.001117 

0.022496 

-0.006590 

0.023754 

-0.017301 

0.023208 

-0.03 1153 

0.019973 

-0.04821  4 

0.012924 

-0.068537 

0.000652 

-0.092191 

-0.01  8494 

-0. 1 19527 

-0.046586 

-0.1 42539 

-0.074273 

-0. 1 19527 

-0.046586 

-0.092191 

-0.018494 

-0.068537 

0.000652 

-0.048214 

?.01  2924 
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-0*031153 

0.019973 

-0.0  17301 

0.  02  3208 

-0*006590 

0.023754 

0.001117 

0.022496 

0.006064 

0.020127 

0.008605 

0.017170 

0.009185 

0.014005 

0.00831 2 

0.010879 

0.006510 

0.007930 

0.004275 

0.005203 

0.002036 

0.00  2664 

AXIAL  MAGNETIC  CURRENT 

ON  THE  SLOT= 

0. 1 79188 

-0.096876 

0.243137 

-0.  133158 

0.325701 

-0. 179775 

0.392747 

-0.21  8353 

0.454081 

-0.253917 

0 .50661 5 

-0.  28  4658 

0.550593 

-0.310557 

0.585321 

-0.331120 

0.610444 

-0.34  6054 

0.625641 

-0.  355113 

0.630727 

-0.358149 

0.625640 

-0.355113 

0.610443 

-0.346054 

0.565320 

-0.331120 

0. 550592 

-0.310557 

0.506615 

-0.284657 

0.454080 

-0.  253917 

0.392747 

-0.218353 

0.325700 

-0.  179775 

0.243136 

-0. 133157 

0.  179188 

-0.096876 

THE  EQUIVALENT  NETWORK 

PARAMETERS 

y=  0.96  1  138E— 02 

0.82201 8E— 02 

V0=  -0.951908E  02 

0.8 18828E  02 

the  power  transmitted  p 

=  0.131187E 

REFERENCES 


[1]  D.  Kajfez,  "Excitation  of  a  Terminated  TEM  Transmission  Line  Through 
a  Small  Aperture,"  AFWL  Interaction  Note  215,  July  1974. 

[2]  C.  M.  Butler  and  K.  R.  Umashankar,  "Electromagnetic  Excitation  of  a 
Wire  Through  an  Aperture-Perforated  Conducting  Screen,"  IEEE  Trans. , 
vol .  AP-24 ,  pp.  456-462,  July  1976. 

[3]  K.  R.  Umashankar  and  J.  R.  Wait,  "Electromagnetic  Coupling  to  an 
Infinite  Cable  Placed  Behind  a  Slot-Perforated  Screen,"  IEEE  Trans. , 
vol.  EMC-20,  No.  3,  pp.  406-411,  August  1978. 

[4]  R.  F.  Harrington,  Time-Harmonic  Electromagnetic  Fields,  McGraw-Hill 
Book  Co.,  New  York,  1961. 

[5]  R.  F.  Harrington,  Field  Computation  by  Moment  Methods,  The  Macmillan 
Co.,  New  York,  1968. 

[6]  R.  F.  Harrington  and  J.  R.  Mautz,  "A  Generalized  Network  Formulation 
for  Aperture  Problems,"  IEEE.  Trans.,  vol.  AP-24,  pp.  870-873, 
November  1976. 

[7]  C.  M.  Butler  and  D.  R.  Wilton,  "General  Analysis  of  Narrow  Strips 
and  Slot,"  IEEE  Trans.,  vol.  AP-28,  pp.  42-48,  January  1980. 

[8]  H.  H.  Chao  and  B.  J.  Strait,  "Computer  Programs  for  Radiation  and 
Scattering  by  Arbitrary  Configurations  of  Bent  Wires,"  Scientific 
Report  No.  7  on  Contract  No.  F19628-68-C-180,  Syracuse  University, 
September  1970. 

[9]  R.  F.  Harrington,  "Matrix  Method  for  Field  Problems,"  Proc.  IEEE, 
vol.  55,  pp.  136-149,  February  1967. 

[10]  J.  Luzwick  and  R.  F.  Harrington,  "Computer  Program  for  Mutual 
Coupling  in  a  Finite  Planar  Rectangular  Waveguide  Antenna  Array," 
Technical  Report  No.  8,  on  Contract  No.  N00014-76-C-0225, 

Syracuse  University,  July  1978. 

[11]  R.  F.  Harrington  and  J.  R.  Mautz,  "Reactively  Loaded  Directive 
Antennas,"  Technical  Report  TR-74-6,  Syracuse  University, 

September  1974. 


